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ABSTRACT 


This  memorandum  records  the  proceedings  of  a  two-part  seminar  series  sponsored  by  the 
Naval  Undersea  Warfare  Center  Division,  Newport,  RI,  during  1993.  The  first  part  of  the 
seminar  series  dealt  with  the  application  of  microfabrication  to  turbulence  control;  the  second 
part  dealt  with  vortex  interactions. 
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FOREWORD 


The  Naval  Undersea  Warfare  Center  Division,  Newport,  RI,  sponsored  a  two-part  seminar 
series  in  1993  titled  Microfabrication  and  Turbulence.  The  first  part  of  the  series  comprised  four 
seminars  devoted  to  the  application  of  microfabrication  to  turbulence  control;  the  second  part 
comprised  three  seminars  that  dealt  with  vortex  interactions.  The  seminars — held  in  February, 
April,  May,  June,  and  July — were  an  activity  of  NUWC  Newport's  Hydrodynamics  Sphere  of 
Excellence ,  which  is  one  of  NUWCs  leadership  areas.  The  presentation  materials  used  during 
the  seminars,  consisting  mainly  of  viewgraphs,  are  reproduced  in  this  report  in  their  original 
form. 


In  the  first  seminar,  Professor  Kenneth  Breuer  of  MIT  used  scaling  arguments  to  show  that 
only  microfabrication  offers  the  technology  to  implement  active  control  strategies  in  a  turbulent 
boundary  layer. 

In  the  second  seminar,  Professor  Gad-el-Hak  of  the  University  of  Notre  Dame  reviewed 
various  turbulence  control  methodologies  and  gave  order-of-magnitude  arguments  to  compare 
the  penalties  and  benefits  inherent  in  each. 

While  many  electrical  engineers  are  developing  the  technology  of  microfabrication, 
Professor  Chih-Ming  Ho  of  UCLA  is  examining  the  science  of  the  subject.  In  the  third  seminar, 
he  compared  the  theory  and  experiments  of  the  fluid  mechanics  of  the  process  of  etching  micron¬ 
sized  channels. 

Professor  Yu-Chong  Tai  of  CalTech  presented  the  fourth  seminar,  in  which  he  reviewed 
microfabrication  technology  and  then  demonstrated  a  spectrum  of  applications  ranging  from 
neuron  farms  and  micron-sized  wind  tunnels  to  the  development  of  a  host  of  fluid  dynamics 
sensors  of  unprecedented  resolution. 

In  the  second  part  of  the  seminar  series,  vortex  interaction  in  turbulence  was  examined. 

The  turbulence  production  process  in  a  high-Reynolds-number  practical  turbulent  boundary  layer 
is  poorly  understood,  and  it  involves  a  diverse  range  of  scales.  To  complicate  matters,  even  at  a 
very  low  Reynolds  number,  a  typical  coherent  structure  dimension  can  vary  by  a  factor  of  as 
much  as  1 0.  To  understand  the  regeneration  process  of  the  constituent  hairpin  vortices,  Professor 
C.  R.  Smith  of  Lehigh  University  has  conducted  very-low-Reynolds-number  experiments  on 
what  can  be  termed  turbulence  with  lesser  variability.  In  the  fifth  seminar,  Professor  Smith 
showed  that  his  approach  has  elucidated  several  features  of  the  production  process,  particularly 
the  emergence  of  a  thin,  vertical-shear  layer  near  the  wall.  He  also  presented  the  theoretical 
works  of  Professor  J.  D.  A.  Walker,  which  complement  his  experiments. 


v 


In  the  sixth  seminar,  Professor  Demetri  Telionis  of  Virginia  PolyTech  discussed  the 
strategies  of  "chopping"  discrete  coherent  vorticities  and  reducing  the  pressure  spikes  caused  by 
the  vorticities  impacting  the  leading  edge  of  a  hydrofoil. 

The  final  seminar  of  the  series  was  given  by  Professor  A.  Leonard  of  CalTech.  In  the 
absence  of  a  wall,  in  two-dimensional  and  axisymmetric  flows,  his  vortex  methods  closely 
simulate  several  complex  flow  experiments.  However,  Professor  Leonard  pointed  out  that, 
because  of  three-dimensional  viscous  effects  calculating  drag  due  to  the  flow  past  a  seemingly 
two-dimensional  flat  plate  in  the  presence  of  a  grazing  wall  remains  a  challenge. 
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Part  I: 


Application  of  Microfabrication  to  Turbulence  Control 


1-1.  Micromachined  Sensors  and  Actuators 
for  the  Measurement  and  Control 
of  Turbulent  Boundary  Layers 

Kenneth  $.  Breuer 

Massachusetts  Institute  of  Technology 
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Micromachined  sensors  and  actuators  for  the 
measurement  and  control  of  turbulent  boundary 

layers 

by 

Kenneth  S.  Breuer 

Department  of  Aeronautics  and  Astronautics 
Massachusetts  Institute  of  Technology 


The  emerging  field  of  micromachining  has  opened  a  new  door  to  a  wide 
variety  of  sensors  and  actuators  that  can  be  designed  to  implement  bound¬ 
ary  layer  control.  Using  manufacturing  techniques  borrowed  from  the  mi¬ 
croelectronics  industry,  one  can  construct  micro-electro-mechanical  systems 
(MEMS)  for  the  measurement  of  pressure,  shear  stress,  etc.  as  well  as  ac¬ 
tuators  to  provide  closed  loop  flow  control.  These  devices  are  small  (50-100 
microns),  robust,  high-bandwidth  and  very  cheap  and  seem  ideal  for  deploy¬ 
ment  in  massive  arrays  for  distributed  flow  sensing  and  control.  Examples  of 
the  state-of-the-art  in  MEMS  for  flow  control  are  presented  along  with  some 
of  the  problems  that  must  still  be  addressed. 

Dynamical  issues  associated  with  active  flow  control  are  also  discussed. 
The  structure  and  dynamics  of  the  “most  dangerous”  elements  in  a  turbulent 
boundary  layer  are  presented  along  with  possible  techniques  for  their  control 
using  massively  distributed  sensor- actuator  arrays.  Some  comments  on  the 
application  of  this  technology  to  the  control  of  transition  to  turbulence  are 
also  offered. 


Thursday,  the  25th  February  1993 
Conference  Room  B,  Bldg.  102 
Time:  10:30  AM 
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Micromachined  sensors  and  actuators  for  boundary 

layer  control 

Kenneth  S.  Breuer 

Department  of  Aeronautics  &  Astronautics 
Massachusetts  Institute  of  Technology 

Seminar  at  NUWC,  Newport,  RI 
March  2, 1993 


(With  assistance  from  Marty  Schmidt  &  Errol  Arkilic) 


Outline 


Review  of  near-wall  turbulent  boundary  layer  structure 
Dynamical  model  for  coherent  structure  evolution  and  control 
Limitations  on  practical  control 
Review  of  Micromachined  sensors  and  actuators 
Application  of  micromachining  for  boundary  layer  control 
Problems  with  micromachined  devices 


Near- wall  turbulent  boundary  layer  structure 

Coherent  Structures 


Large  fraction  of  rms  wall-pressure  fluctuations  associated  with  large-amplitude  pressure  peaks 

"Waves"  observed  (?)  in  pressure  and  vertical  velocity 
-  instablity/breakdown  of  coherent  structure  (?) 


Near-wall  turbulent  structure 

Inclined;  shear  layers 


Reynolds  stress  and  pressure  peaks  associated  with  elongated,  inclined,  shear  layer  structures 
Non-symmetric  structures  predominate  (a  =  0  is  present) 


Dynamical  model  for  near  wall  turbulent  structures 
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•  Main  evolution  dominated  by  turbulent  mean  shear  flow 

•  Viscous  and  nonlinear  effects  take  place  on  a  much  longer  (unimportant)  timescale 


MO 


Dynamical  model  for  near  wall  flows 

a  =  0  component 
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•  Vertical  component  can  (and  does)  decay 

•  Streamwise  component  grows  algebraically  due  to  three-dimensional  "liftup"  effect 


Numerical  simulation  of  coherent  structures 

Initial  conditions 


Two  pairs  of  counter-rotating  vortides 
Result  of  a  localized  burst  of  Reynolds  stress 
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Numerical  simulation  of  coherent  structure  evolution 

(linear  inviscid  model) 


.x/X* 


•  Development  of  inclined  shear  layer 

•  Intensification  as  time  advances 

•  (No  saturation  or  breakdown  due  to  model  simplicity  -  not  important  for  control) 


M3 


Methods  for  turbulent  flow  control 


U(y) 


•  Model  active  wall  as  localized  wall  motion  (all  wavenumbers  excited) 

•  Dominant  effects  are: 

-  a  =  0;  dominant  algebraic  growth  term 

-  Lag  due  to  non-collocation  of  sensor  and  actuator 

-  Assume  spanwise  motion  decoupled  with  z+  ~  100 
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Problems  with  the  implimentation  of  turbulent 

boundary  layer  control 

(Time  scales) 


Minimum  Turbulent  Time  Scales 


Unit-Re  (1/ft) 
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Problems  with  the  implimentation  of  turbulent 

boundary  layer  control 

(Length  scales) 


Minimum  Turbulent  Length  Scale 


Unit  Re  (1/ft) 
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Development  of  micro-sensors  and  actuators 

Example:  Micro  Hot-wire 


Alignment  pins 
(for  stacking) 


On-chip  electronics 
SCALE:  lOOum 


Hot  film 


•  Sensing  element  size  0.5  x  100  pm 

-  Thin  film  polysilicon  bridge  (good  mechanical  strength,  selectable  resistance) 

•  Stackable  design  -  form  multi-probe  "rakes"  for  boundary  layer  measurements 

-  10  probes  in  a  0.5  mm  span 

-  Necessary  for  pressing  need  for  high-Reynolds  number  instrumentation 


1-17 


What  is  Micromachining 

What  are  its  advantages 


•  Use  of  VLSI  circuit  manufacturing  techniques  for  fabrication  of 
Micro-electromechanical  (MEM)  devices: 

—  Mechanical  devices:  gears,  levers,  micro-honeycombs,  etc 

—  Sensors:  pressure,  shear  stress,  temperature,  chemical,  etc 

—  Actuators:  pressure,  force,  strain,  motors,  etc 

•  Small  dimension:  0(lfim)  in  horizontal  direction,  O(IOOA)  in  vertical 
direction. 

•  Ability  to  integrate  on-chip  electronics  -  reduce  noise,  lab  complexity. 

•  High  Frequency  response:  DC  -  500  kHz,  limited  by  sensitivity  and  electronics. 

•  Cheap:  once  commercialized,  can  be  manufactured  by  the  thousand  at  cost  of 
“pennies”  per  device  (e.g.  DELCO  pressure  sensor). 
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Fabrication  Techniques 


A  brief  overview 


•  Start  with  single  crystal  silicon  wafer 

•  Build  devices  by  cutting  into  wafer  (“Bulk  micromachining),  putting  materials 
on  the  surface  ( “Surface  micromachining)  and  gluing  wafers  together  (Wafer 
bonding) 

•  Use  “Masks”  to  selectively  protect  and  expose  various  layers  to  each  operation 

•  Replicate  the  pattern  across  the  wafer  (10,000  elements  on  a  single  4-inch 
wafer) 

•  Cut  up  the  wafer 

•  Package  (not  a  trivial  step) 

•  Use/S  ell. 
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Fabrication  Techniques 

(continued) 


•  Additive  techniques 

—  Epitaxy  -  crystal  growth 


v  x  v  S  r  N  \  ^  \  \  \  \  v  E>'  P  •  ••■) 

///  /  /  /  /  /  z7  /  /  S,u,CON 

Chemical  Vapor  Deposition  (CVD)  -  addition  of  thin  layers  of  almost  any 
material  ( Si02 >  SizN$,  Poly -Si,  metals,  etc). 


b  * 


. _  .w 


Doping  -  Implantation  of  p+,  n+  silicon,  Boron,  Phosphorus,  etc  for: 
(i)  Electronic  circuits,  (ii)  Etch-stop  control 

\\>\>'  W\\  \  \  ^  S^uaC^M 
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Fabrication  Techniques 

(continued) 


•  Subtractive  techniques 

—  Isotropic  etching  of  Silicon 


js: — 


s; 


*0  VJ  (y ! 

(^VVW  sS 


No  spatial  dependence,  typically  a  HNOz,HF  bath,  masked  by  SizN±  (or 

Si02). 


Fabrication  Techniques 

(continued) 


Subtractive  techniques  (continued) 
—  Anisotropic  etching  of  Silicon. 


Follows  Crystal  boundaries.  Very  precise  control  of  shape  and  geometry. 

KOH  bath,  can  be  masked  by  Si02 

Can  be  used  to  machine  bridges  and  cantilevers. 
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Fabrication  Techniques 

(continued) 


•  Subtractive  techniques  (continued) 

—  Isotropic  etching  of  other  materials  (Sz02,  Si^N^  Poly -Si,  metals,  etc)  - 
All  kinds  of  etches  and  resists  available. 

•  Wafer  Bonding 

—  Silicon  equivalent  of  welding  -  Extremely  strong  joint  formed. 

—  in  conjunction  with  doping  techniques,  can  produce  very  accurate 
dimensional  control  in  the  vertical  direction. 


\ 

‘  CAvrry 
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Examples  -  Pressure  Sensor 

Readout,  options 


Simple  and  accurate  but  has  temperature  dependence 
•  Capacitive  readout: 


More  complex,  requires  on-chip  electronics,  but  very  sensitive. 
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Microfabricated  Floating-Element  Schematic 


Si  Epi:  Layer 


Piezoresist've  Tethers 


Piezoresistive  Readout  Scheme 


Flow  Direction 


Si  Epi  Layer 


•  Linear  Response 

•  Dielectric  Isolation 

=>  High  Temperature  Operation 


AVout 
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p++  Region 


Summary  of  the  Shear-Stress  Sensor 


Top  View 


f 
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Fabrication  Summary  (cont.) 


Si  Epi  Layer 


Si  Epi  Layer 
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Fabrication  Summary  (cont.) 

_ Top  View _ 


PECVD  Oxide 
Amorphous  Si 
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Fabrication  Summary  (cont.) 


Top  View 
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Micromachined  Actuators 

Preliminary  designs 


•  Membranes  (high  forces,  low  displacements) 

•  Resonant  cantelevers  (low  forces,  high  displacements) 
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Problems  with  microsensors  and  actuators 

Micro-fluid  mechanics 


•  Fluid  mechanics  at  micron  scale  is  not  completely  straightforward 

•  Gases:  Mean-free-path  effects  become  significant  (High  Knudsen  number,  low  Reynolds 
number) 

-  Viscous,  transonic  flow! 

-  anomolously  low  viscosities 

•  Liquids:  double  layers  and  dielelctric  effects  (especially  with  polar  liquids) 

-  more  anomolous  viscosities,  charges,  etc 
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Conclusions 


•  Wall  pressure  fluctuations  dominated  by  passage  of  near-wall  turbulent  structures 

•  Near-wall  coherent  structures  characterized  by  inclined  shear  layers. 

-  dominated  by  long-wave  structures  (a  =  0) 

•  Dynamics  dominated  by  "algebraic  instability"  (linear,  inviscid,  three-dimensional) 

•  Simple  theory  leads  to  dynamical  equations  for  evolution  and  control 

•  Scales  of  normal  turbulent  flows  require  micro-fabricated  sensors  and  actuators 

•  Micromachining  leads  to  promise  of  "smart  skins"  for  boundary  layer  (and  transition)  control 

•  Questions  still  remaining: 

-  Micro-fluid  mechanics  (gaseous  and  liquid  behavour  at  micron  scale) 

-  Signal  processing  for  sensor/ actuator  arrays 

-  Control  algorithms  for  massively  distributed  systems 


1-2.  Interactive  Control  of  Turbulent  Boundary  Layers: 

A  Futuristic  View 


Mohamed  Gad-el-Hak 
University  of  Notre  Dame 
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INTERACTIVE  CONTROL  OF  TURBULENT  BOUNDARY  LAYERS: 

A  FUTURISTIC  VIEW 

Mohamed  Gad-el-Hak 

Department  of  Aerospace  &  Mechanical  Engineering 
University  of  Notre  Dame 
Notre  Dame,  IN  46556 


ABSTRACT 

The  ability  to  actively  or  passively  manipulate  a  flow  field  to  effect  a  desired 
change  is  of  immense  technological  importance.  In  its  broadest  sense,  the  art  of 
flow  control  probably  has  its  roots  in  prehistoric  times  when  streamlined  spears, 
sickle-shaped  boomerangs,  and  fin-stabilized  arrows  evolved  empirically  by 
archaic  Homo  sapiens.  In  this  presentation,  methods  of  control  to  achieve 
transition  delay,  separation  postponement,  lift  enhancement,  drag  reduction, 
turbulence  augmentation,  or  noise  suppression  are  considered.  An  attempt  is 
made  to  present  a  unified  view  of  the  means  by  which  different  methods  of 
control  achieve  a  variety  of  end  results.'  It  is  argued  that  future  systems  for 
control  of  turbulent  flows  in  general  and  turbulent  boundary  layers  in  particular 
would  greatly  benefit  from  the  merging  of  a  new  science  and  a  novel  technology. 
Control  of  chaotic,  nonlinear  dynamical  systems  has  been  demonstrated 
theoretically  as  well  as  experimentally,  at  least  for  simple,  one-degree-of- 
freedom  systems.  Microfabrication  is  an  emerging  technology  which  has  the 
potential  for  producing  inexpensive,  programmable  sensor/actuator  chips  that 
have  dimensions  of  the  order  of  a  few  microns.  Together,  these  two  disciplines 
could  provide  efficient  interactive  control  systems  to  achieve  a  desired  beneficial 
goal  for  practical  flow  fields. 

Thursday ,  the  8th  April  1993 
Conference  Room  B,  Bldg.  102 
Time:  10:30  AM 
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INTERACTIVE  CONTROL  OF  TURBULENT  BOUNDARY  LAYERS: 

A  FUTURISTIC  VIEW 


Mohamed  Gad-el-Hak 

Department  of  Aerospace  &  Mechanical  Engineering 
University  of  Notre  Dame 
Notre  Dame,  IN  46556 

ABSTRACT 

The  ability  to  actively  or  passively  manipulate  a  flow  field  to  effect  a 
desired  change  is  of  immense  technological  importance.  In  its  broadest  sense,  the 
art  of  flow  control  probably  has  its  roots  in  prehistoric  times  when  streamlined 
spears,  sickle-shaped  boomerangs,  and  fin-stabilized  arrows  evolved  empirically 
by  archaic  Homo  sapiens.  In  this  presentation,  methods  of  control  to  achieve 
transition  delay,  separation  postponement,  lift  enhancement,  drag  reduction, 
turbulence  augmentation,  or  noise  suppression  are  considered.  An  attempt  is 
made  to  present  a  unified  view  of  the  means  by  which  different  methods  of 
control  achieve  a  variety  of  end  results.  It  is  argued  that  future  systems  for 
control  of  turbulent  flows  in  general  and  turbulent  boundary  layers  in  particular 
would  greatly  benefit  from  the  merging  of  a  new  science  and  a  novel  technology. 
Control  of  chaotic,  nonlinear  dynamical  systems  has  been  demonstrated 
theoretically  as  well  as  experimentally,  at  least  for  simple,  one-degree-of- 
freedom  systems.  Microfabrication  is  an  emerging  technology  which  has  the 
potential  for  producing  inexpensive,  programmable  sensor/actuator  chips  that 
have  dimensions  of  the  order  of  a  few  microns.  Together,  these  two  disciplines 
could  provide  efficient  interactive  control  systems  to  achieve  a  desired  beneficial 
goal  for  practical  flow  fields. 
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Aim 


Flow  control: 

The  ability  to  actively  or  passively 
manipulate  a  flow  field  in  order  to 
effect  a  beneficial  change. 


Provide  a  broad  overview  of  the 
arts  and  sciences  of  flow  control. 


Focus  on  control  of  wall-bounded 
flows: 

Immense  technological  importance. 
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Aim 


Control  of  a  TBL  to  achieve 
a  variety  of  beneficial  changes. 


Unifying  principles? 


Coherent  structures. 


Selective  control. 


Outlook  for  the  future. 
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Past 


Archaic  Homo  sapiens: 

Streamlined  spears;  sickle-shaped  boomerangs; 
fin-stabilized  arrows. 


Present: 


Prandtl  (1904)  pioneered  the  science  of  flow  control. 

Boundary  layer  theory  has  its  roots  in  a  question  that  boggled 
Prandtl  for  few  years: 

Failure  of  conical  diffusers  in  a  large  air  duct 
to  achieve  the  expected  pressure  recovery. 


<*“  Future: 


Empirical  Art 


Predictive  Science 


age  CFD 


%  Nonlinear  Dynamical 
Systems  Theory 

+ 

%  Microfabrication 


Interactive  control 
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15  x  109 
years  ago 


Big  Bang 


4.6  x  109 


Solar  System  formed 


1.5  x  109 
0.25  x  109 
4  x  106 


Life  began  on  Earth 
Dragonfly 

Ancestral  hominid 


2  x  106 


Homo  habilis  (handy) 


1.5  x  10* 

500,000 

200,000 

10,000 

years  ago 


Homo  erectus 

Archaic  Homo  sapiens 

Homo  sapiens 

Agriculture 

Present 


(wise) 
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10,000 

years  ago 


Agriculture 


5.500 

5,000 

3.500 

300 

150 

<  100 
years  ago 

Present 


Plough 


Irrigation  Systems/Writing 


Metals 


Modern  Science/Newton's  Laws 


Navier-Stokes  Equations 


Powered  Flight/ 
Boundary  Layer  Theory/ 
Science  of  Flow  Control 


Future 


97999777 


Flow  Control  Goals 

Transition  delay/advancement. 

Turbulence  enhancement/ 
suppression/relaminarization. 

*•*  Separation  prevention/provocation. 

II  II  II  II 

*•"  Skin-friction/pressure  drag  reduction. 
Lift  enhancement. 


«•*  Heat  transfer/mixing/ 

chemical  reaction  augmentation. 


«•“  Noise  suppression. 
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Desired  Traits  for  a 
Flow  Control  Device 


Simple. 

Inexpensive  to  build. 
Inexpensive  to  operate. 
Has  no  trade-offs. 

^  Does  not  exist! 
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Drag  Reduction 


Pressure  Drag: 

*  Separation 

*  Displacement 

*  Wave 

*  Induced 


Skin  Friction: 

*  LBL 

*  Transition 

*  TBL 
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Drag  Reduction 


Internal  Flows:  Close  to  100%  is  skin  friction. 
^  Throughput;  pumping  power;  duct  size/cost. 


External  Flows:  Land,  air,  and  sea  vehicles. 
♦  Range;  fuel  cost/volume;  payload;  speed. 


) 


Submarines:  90%  skin  friction. 

Aircraft:  50%  skin  friction;  25%  induced  drag. 

Annual  fuel  costs  for  U.S.  commercial 
fleet  exceed  $10B. 

It  costs  $2B  to  develop  a  new  wing. 
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Tools  for  Controlling 


Surface: 

*  Roughness;  Riblets;  Fences 

*  Curvature 

*  Shape 

*  Compliant 

*  Mass  transfer 

*  Heat  transfer 


Freestream: 

*  LEBU 

*  Acoustics 

*  Turbulence  levels;  Gust 

Additives 

*  Polymers;  surfactants 

*  Micro-bubbles 

*  Particles;  dust;  fibers 
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Types  of  Control  Devices 


Passive; 

■  No  energy  expenditure 

■  e.g.,  riblets 


Active: 

■  Requires  energy 

■  e.g.,  suction 

■  Is  saving  larger  than  penalty ? 

^  Interactive: 

©  Control  input  based  on  measurement 
@  Open  or  feedback  control  loop 
©  Sensor(s)  +  actuator(s)  needed 
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Is  there  a  unifying  principle 
that  ties  together  a  majority 
of  flow  control  devices? 


Help  to  understand  the  physics. 


Help  in  contemplating  new  artifacts. 


Help  in  minimizing  adverse  tradeoffs. 
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Vorticity  Framework 


Wall  as  a  source/sink  of 
spanwise/streamwise  vorticity. 

Magnitude  and  direction  of  vorticity  flux. 


u. 


*  Lighthill  (1963): 
as  a 

*  Willmarth  (1975): 

*  Reshotko  (1976): 

&  Reynolds&Carr  (1985): 


Place  boundary  layer 

correctly  in  the  flow 
whole. 

Nature  and  structure  of 
TBL. 

Transition  control. 
Separation  control. 
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1-54 


-56 


U(y) 


■  Suction 

■  Favorable  P-grad. 

■  Heating  (water) 


U(y) 
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Velocity  Profile  Modifiers 


A  fuller  velocity  profile  (negative 
curvature  or  vorticity  source  at  wall) 
will  lead  to: 

separation/ transition  resistance; 

turbulence  suppression; 
higher  skin  friction. 


*  Wall  suction;  vw  <  0. 


*  Favorable  pressure 
gradient;  dP0/dx  <  0. 

*  Lower  viscosity  at  wall; 
9ji./3y  >  0. 
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LBL:  Cq  ~  Re'0-5 

To  prevent 
separation 

TBL:  Cq  =  0.002  -  0.004 


dSe 

(  Cf  =  2  +  2  Cq 

+  Pressure  gradient  terms  ). 
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At  Re  =  107:  Cf|T  =  0.003 

Cf|L  =  0.0004 

-  Cf  =  0.0030 


Suction 


Cf  =  0.0004 
Cf  =  0.0000 


<•*  To  delay  transition: 

Laminar  separation: 

<•"  Turbulent  separation: 

Asymptotic  profile  (TBL): 

**  Relaminarization: 


Cq  =  0.0001 
Cq  =  0.0003 
Cq  =  0.003 
Cq  =  0.003 

Cq  =  0.01 
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Passive  Suction 


I  In  high-speed  flows. 

I  Porous  surface  underneath  shock 
wave. 

I  Mass  self-bleed  from  downstream 
of  shock  to  upstream  to  mitigate 
the  local  pressure  gradient. 


u 


Shock  Wave 


i 

i  i 

L^/J 

n 

f 

A  series  of  weaker  shock  waves; 

prevent  separation; 
reduce  wave  drag. 
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Shaping 


LBL:  If  fluid  decelerates  faster  than 

Uo  -  x-  °-09 

boundary  layer  separates. 


TBL:  If  fluid  decelerates  faster  than 

U0  -  x-  °-23 

boundary  layer  separates. 


Efficient  momentum  transport  provides 
the  mechanism  for  mixing  slower  fluid 
near  wall  with  faster  fluid  particles 
further  out. 


Ife  Transition  advance. 
Turbulators. 
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Stratford  (1959)  Closure 


Flow  with  continuously  zero  skin- 
friction  throughout  its  region  of 
pressure  rise. 


Achieve  a  specified  pressure  rise  in 
the  shortest  possible  distance  and 
with  the  least  possible  dissipation 
of  energy. 


Apply  Stratford  closure  immediately 
after  transition  to  maximize  L/D. 


Liebeck  (1978) 

Assume  incipient  separation 

Calculate  pressure  field  required 

Use  inverse  calculation 
procedure  to  derive  shape 

Major  Concern: 
Off-design  conditions 
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Hybrid  Laminar  Flow 

Control 

Suction  at  leading  edge  region 

+  careful  contouring  of  airfoil 
shape  elsewhere. 

Very  promising  to  delay  separation 
as  well  as  to  prevent  transition. 


High  lift  during  takeoff  and  landing. 
Lower  skin  friction  during  cruise. 
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Lower  viscosity  at  wall 


Cooling  wall  in  gases. 

(also  increases  density  and,  thus,  near-wall 
momentum) 

^  Heating  wall  in  liquids. 

Film  boiling. 

Cavitation. 

Sublimation. 

Chemical  reaction. 

*•“  Wall  injection  of  a  lower-viscosity 
fluid. 

Shear-thinning  additives. 

) 
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Wall  Heating/Cooling 


Lower  wall  viscosity  leads  to: 

separation/transition  resistance; 

turbulence  suppression; 
higher  skin  friction. 


Demonstrated  for  high-speed  flows 
in  gases 


♦  Requires  a  heat  sink 

(e.g.  Cryognically-fueled  aircraft) 


In  liquids,  wall  heating  lowers  near 
wall  viscosity  but  density  remains 
essentially  unchanged. 

#  Very  little  research 

^  How  about  marine 
propellers? 


New  Way  of  Looking  at 
Turbulent  Flows 


Coherent  structures  seem  to  be 
ubiquitous  in  all  turbulent  shear 
flows. 


**  Enormous  technological  interest 
in  that  they  offer  the  possibility  for 
controlling: 

Jets;  wakes;  mixing  layers; 
wall-bounded  flows. 


In  a  TBL,  hierarchy  of  quasi- 
periodic  structures: 

Large  eddies;  Falco’s  eddies; 
low-speed  streaks;  bursts. 
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Little  use  of  structural 
information  in  searching  for 
innovative  control  methods. 


Possible  reason: 


Flow  control  is  as  old  as  man, 
while  organized  structures  were 
shown  to  exist  in  a  TBL  only 
30  years  ago. 


At  best,  we  are  presently  trying  to  use 
structural  information  to  explain  why  a 
specific  control  device  works. 


^  How  about  trying  to  create 
control  artifacts  to 
suppress/enhance 
a  particular  structure? 
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COHERENT  STRUCTURES 


Large  Outer  Structures. 


Intermediate  Falco  Eddies. 


Near-Wall  Events. 


&  Low-speed  streaks. 
&  Ejection. 

&  Sweep. 


Bursting 
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Polymers;  Particles;  Riblets 


#  Increased  losses  in 
near-wall  region. 

O  To  stay  in  equilibrium: 
Energy  gain  goes  up. 

C>  Scales  go  up: 

gain  ~  (size)2 
loss  ~  size 

O  Buffer  layer  thickens. 
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C>  Position  for  maximum 

u2  and  uv  farther 
away  from  wall. 


O  Production  of  Reynolds 
stress  inhibited. 


O  Momentum  transport 
lowered. 


O  Skin  friction  lowered. 
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Is  skin-friction  reduction 
associated  with  turbulence 
suppression? 


($5  X  c!!)« 

Polymers;  Particles;  LEBUs; 
Riblets. 


[Act  selectively  on  a  particular  structure] 


&  $5  No: 


Suction;  Wall  Cooling/Heating; 
Favorable  Pressure  Gradient. 


[Act  globally  on  all  eddies] 
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Suction,  etc.,  acts  directly  on 
mean  flow. 


Mean  velocity  modifiers. 
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SUCTION 


■  Flat  Plate: 

Cf  =  2(d50/dx)  +  2Cq 

U  Ji  U 

II  JJ,  u 

&  No  suction: 

0.003  2x0.015  0.0 

Suction  (asymptotic  velocity  profile)! 

0.006  0.0  2x0.003 
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Control  of  a  TBL 


Global. 

Selective: 

#  By  the  flow. 


#  By  design. 


Near -wall  events: 


Very  intermittent  and  random 
in  space  &  time. 

Temporal  phasing  and 
spatial  selectivity  are 
needed  for 
selective  control. 
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WHAT  TO  TARGET? 


Low-speed  streaks 


are  the  most 

visible 

reliable 

detectable 

indicators 


of  the  preburst 
turbulence  production 
process. 
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EXAMPLE 


Submarine 


&  U0  =10  m/s 


#  Re  =  107/  m 

#  v/u*  =  2.6  (I 


Aircraft 

(10  km) 


300  m/s 
107/  m 
2.6  |i 
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SENSORS/ACTUATORS 


#  100  wall  units  spanwise. 
(260  p) 

&  1000  wall  units  streamwise. 
(2.6  mm) 


■  Number  of  elements  =  1.5  x  106/m2. 


■  Frequency  =  600  Hz  (submarine). 

=  18  kHz  (aircraft). 


Actuator's  Response: 


©  Wall  displacement  =  26  (I. 

O  Cq  =  0.0006  (Cf  =  0  +  2  x  0.0006) 

©  AT  =  2°C  (heating;  water) 

=  40°C  (cooling;  air) 
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Vision  for  a  Control  System 


Checkerboard  of  wall  sensors  +  actuators. 


♦  Sensors: 

Pressure;  velocity;  wall  shear;  etc. 


^  Actuators: 

Heating/cooling;  suction/injection; 
wall  movement;  etc. 


For  example: 

©  Piezoelectric  devices  under  flexible  skin. 
©  Terfenol-d  materials. 


John  Lumley  (1991). 


Can  it  be  done? 


Breakthrough  #  1: 


Microfabrication 


Breakthrough  #  2: 


Control  of  Chaos 
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Microfabrication 


Micro- Machinery:  | 

Manufacturing  processes  that  can  create 
microscopic  machinery. 


Combine  electronic  and  mechanical 
components  on  a  single  chip: 

®  Sensor/Actuator. 


Sensors: 

@  Pressure;  temperature;  velocity;  mass  flow; 
sound;  smell. 


Actuators: 

@  Motors;  pumps;  moving  rods; 
valves;  gears;  tweezers. 


Micro-robot  on  a  single  chip: 

©  Sensing;  guidance;  and  control  systems. 
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Microfabrication 


Combine  electronic  and 

mechanical  components  on  a 
single  chip: 

Sensor/ Actuator. 


I  Small  size. 

I  Smart  machine. 
I  Inexpensive. 
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SENSORS 


Sensors  detect  instantaneous  wall 
pressure  or  wall  shear-stress. 


^  Location. 

^  Amplitude. 
*  Phase. 


Example: 

^  Laser-powered  picophones  machined  in 
silicon  using  integrated  circuit  fabrication 
technology. 


O  MIT. 

®  Diameter:  0.2  mm  (200  (i). 

(77  wall  units) 
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ACTUATORS 


Suction  slots 


Wall  displacement 


Heating 


(electric-resistance  heating) 


Cooling 


(thermoelectric  devices 
based  upon  the  Peltier  effect  ) 
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ACTUATORS 


Smart  materials  that  might  allow  a 
wing  to  react  and  adjust  to 
fluctuating  wind  conditions. 

(Grumman  Corporation.) 


Terfenol-d: 


Metal  composite  that  changes  its 
length  when  exposed  to  a 
magnetic  field. 


Nonlinear  material. 


®  (classical  control  impractical.) 

®  (potential  candidate  for  chaotic  control.) 


Able  to  exert  large  forces. 


Example:  Apply  magnetic  field  to  change  the 

length  of  Terfenol-d  inside  a  wing. 
(Change  the  wing's  shape.) 
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ADVANCES  IN  MATERIALS 


Photochromic  glass: 

0  Darkens  in  strong  light. 

Electrochromic  films: 

@  Change  color  when  subjected  to  electric 
current. 

Encapsulated  enzymes  inside  glass: 

0  Change  color  in  the  presence  of  certain 
chemicals. 

Electro-rheological  fluids: 

0  Change  from  liquid  to  near-solid  when 
subjected  to  current. 
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ADVANCES  IN  MATERIALS 


Nitinol  (1962): 

®  A  titanium  and  nickel  alloy  with  memory. 
®  Changes  shape  when  heated/cooled. 

@  Used  in  satellite  antennas. 


Magnetoelastic  ribbons: 

®  Nonlinear  materials  that  change  their 
stiffness  in  the  presence  of  varying 
magnetic  Helds. 


Terfenol-d: 


@  Metal  composite  that  changes  its  length 
when  exposed  to  a  magnetic  Held. 


Fiber-optic  skin  to  measure 
deformation  field  with  a  resolution 

of  0.03  (I. 
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Possible  Pitfalls 


Checkerboard  of 


wall  sensors 

+ 

actuators. 


) 

♦  Information  sensed  incomplete. 

II 

O  Might  be  misinterpreted. 


♦  Checkerboard  actuators  might  be 
less  effective. 


1-89 


Nonlinear  Dynamical 
Systems  Theory 


Butterfly  Effect: 


I  Sensitive  dependence  of 
nonlinear  differential  equations 
on  initial  conditions. 


I  Small  disturbances  can 

radically  alter  a  deterministic 
system  to  yield  rich, 
unpredictable  behavior. 
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Chaos 


«•“  Douglas  Hofstadter: 

"It  turns  out  that  an  eerie  type  of 
chaos  can  lurk  just  behind  a  faqade 
of  order— 

and  yet ,  deep  inside  the  chaos  lurks 
an  even  eerier  type  of  order." 


Chaotic  system  continually  shifts 
from  one  periodic  pattern  to  another, 

creating  the  appearance  of  randomness. 

#  Apparent  randomness  of  a  chaotic 
system  is  only  skin  deep . 


Beneath  the  apparent  unpredictability 
lies  an  intricate  but  highly  ordered  structure. 

)  A  complicated  web  of  interwoven 

patterns  of  regular  (periodic)  motion. 
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Control  of  Chaos 


A  chaotic  system, 
though  apparently  random, 
actually  consists  of  a  number  of 
different  periodic  motions  (orbits). 


System  will  shift  from  one  motion  to 
the  other,  abruptly  and  arbitrarily, 
ad  infinitum. 


O  Reverse  Chaos'. 

♦  Pick  one  of  these  orbits  and  lock 
onto  it. 

Apply  a  tiny  push  at  just  the  right 
moment. 
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Control  of  Chaos 

- 

Like  balancing  a  ball  on  a  saddle. 

o 


^  Saddle  stays  still 

li 

Ball  rolls  off. 

♦  Continuously  adjusted  slightly 
Ball  stays  in  place. 
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Control  of  Chaos 


Theory  for  continuous  systems: 
Ott,  Grebogi  and  York  (1990). 

(University  of  Maryland.) 


Control  the  chaotic  behavior  of  a 
parametrically  driven,  gravitionally 
buckled,  amorphous  magnetoelastic 
ribbon. 

Ditto,  Rauseo  and  Spano  (1990). 

(Naval  Surface  Warfare  Center,  Silver  Spring,  MD.) 


Laminarize  a  chaotic  thermal 
convection  loop  (thermosyphon). 

Singer,  Wang  and  Bau  (1991). 

(University  of  Pennsylvania.) 
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Control  of  Chaos 


Chaotic  attractor  typically  has 
embedded  within  it  an  infinite 
number  of  unstable  periodic  orbits. 


) 

Apply  small  time-dependent 
parameter  perturbations  so  as  to 
stabilize  anyone  of  the  existing 

Orbits  (which  reside  on  strange  attractors). 

♦  Lock  system  into  one  particular 
type  of  repeating  motion. 
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Advantages 


Chaotic  motion 

Periodic  motion. 


A  priori  analytical  knowledge  of  the 
system  dynamics  is  not  needed. 


Could  apply  to  any  chaotic  system 
for  which  a  faithful  Poincare  section 
can  be  constructed. 


^  Effective  even  in  the  presence  of 
small-amplitude  noise  or  small 
modelling  errors. 
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Direct  Trajectories  to  Targets 


Exploit  the  exponential  sensitivity 
of  a  chaotic  system  to  minute 
perturbations. 


Direct  the  system  to  a  desired 
accessible  state  in  a  short  time. 


Effective  even  in  the  presence  of 
small-amplitude  noise  or  small 
modelling  errors. 


) 


Substantially  alter  the  future 
state  by  a  minute  perturbation. 
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Procedure 


Direct  Trajectories  to  Targets: 


Sense  present  state. 


Apply  a  judiciously  chosen 
perturbation. 


System  is  rapidly  directed  towards 
the  desired  state. 
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Turbulent  Boundary  Layers 


Recent  advances  in  dynamical 
systems  theory  provide  opportunities 
for  achieving  efficient  control  of 


TBL. 


By  making  small  time-dependent 
adjustments  to  one  of  the  parameters 
governing  the  system's  behavior. 
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Turbulent  Boundary  Layers 


Need  to  establish  a  conceptual 
framework  for  adapting  York’s 
theory  to  control  (intermittent) 
bursting  events  in  TBL. 


Quasi-periodic  events  modelled  as  a 
dynamical  system. 


Control  about  anyone  of  its  many 
unstable  fixed  points 
(periodic  orbits)  embedded  in  the 
chaotic  attractor. 


Aim  is  to  control  the  bursting  events 
superposed  on  the  fixed  point 
solution  of  the  6-mode  (or  higher) 
dynamical  system  that  describes  the 
near-wall  activities. 
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Bursting  Phenomenon 


Use  proper  orthogonal 
decomposition  to  extract  a 
low-dimensional  dynamical  system 
from  experimental  data  for  the 
near-wall  region  of  TBL. 

Aubry,  Holmes  and  Lumley  (1988;  1991). 
(Cornell  University.) 


**  Expand  instantaneous  velocity  field 
using  experimentally  determined 
eigenfunctions  which  are  in  the  form 
of  streamwise  rolls. 

Apply  Galerkin  projection. 

Truncate  to  obtain  a 
low-dimensional  set  of  (six)  ODE. 

Examine  fixed  points  and  topology  of  their 
stable  as  well  as  unstable  manifolds. 


Chaos  as  well  as 

regular  intermittency  are  exhibited. 
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Control  of  Bursting 


Proper  orthogonal  decomposition. 

(an  optimally  convergent  representation.) 


Dynamical  system  leaves  one  fixed 
point  and  jumps  to  another  along  a 
hetroclinic  cycle. 


Holding  system  near  first  fixed  point 
^Delaying  jump 
-►Drag  reduction. 
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Procedure 


) 


Construct  a  Poincare  map  to  classify 
the  type  of  intermittency 

Xj+i  =  f(Xi,p) 


Xi  discrete  time  signal 
p  parameter  to  be  controlled 


*•*  Apply  a  small  time-dependent 
perturbation  to  the  parameter 
controlling  the  stability  of  the  fixed 
point. 


Force  the  eigenvalue  to  remain 
within  the  circle  of  stability  that 
surrounds  the  (unstable)  fixed  point. 
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Possible  Pitfalls 


System  with  infinite  number  of  degrees  of 
freedom  are  not  readily  susceptible  to  an  easy 
dynamical  systems  approximation. 


Noise  in  the  system  tends  to  kick  the 
orbit  out  of  the  circle  of  stability 

(surrounds  the  unstable  fixed  point). 

^  Forces  the  operator  to  increase  the  control 

amplitude  in  order  to  keep  the  orbit  close  to  the 
fixed  point. 


Manifold  along  which  the  system 
leaves  fixed  point  might  not  be 
one-dimensional. 


A  burst  is  assumed  to  leave  a  fixed 
point  along  the  average  path. 

Actuator  pushes  back  along  same  path. 

In  reality,  most  bursts  would  leave  to  one  side 
or  the  other  of  the  average  path. 
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Microfabrication 


Nonlinear  Dynamical 
System  Theory 


Interactive  control 


Outlook 


Tremendous  energy  saving  potential 
for  vehicles  which  have  notoriously 
high  drag:  Automobiles;  helicopters; 
tractor-trailer  trucks. 


Stand-by  techniques  for  off-design 
situations?? 


Combination  of  approaches?? 

Microfabrication 

+ 

Nonlinear  Dynamical  Systems  Theory 

u 

Interactive  control 


1-106 


Additional  Reading 


Gad-el-Hak,  M.  (1989)  ’’Flow 

Control,”  Applied  Mechanics 
Reviews  42,  pp.  261-293. 


Gad-el-Hak,  M.  (1990)  "Control  of 

Low-Speed  Airfoil  Aerodynamics,” 
AIAA  Journal  28,  pp.  1537-1552. 


Gad-el-Hak,  M.,  and  Bushnell, 

D.M.  (1991)  "Separation  Control: 
Review,”  Journal  of  Fluids 
Engineering  113.  pp.  5-30. 
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Reverse  Blank 


1-3.  Flows  in  Micro  Systems 


Chih-Ming  Ho 

University  of  California  at  Los  Angeles 
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Reverse  Blank 


NUWC  Division  Newport,  R.  I. 
SEMINAR  NOTICE 


FLOWS  IN  MICRO-SYSTEMS 


Professor  Chih-Ming  Ho 
University  of  California,  Los  Angeles 


The  pressure  distribution  of  an  integrated  pressure-sensor/microchannel  sys¬ 
tem  will  be  presented.  We  will  first  discuss  the  scientific  aspects  of  the  etching 
technique  that  is  used  to  manufacture  the  microsystem.  Based  upon  the  findings 
of  the  diffusion  and  the  chemical  reaction  processes,  a  universal  etching  model 
has  been  established.  This  model  can  accurately  predict  the  etching  length  vs. 
time  over  a  wide  range  of  HF  concentrations  (3-49  wt%).  We  then  examine  the 
pressure  distributions  in  the  channels  measured  by  the  integrated  micro  pressure 
transducers.  The  micro  channels  are  one  micron  thick  and  the  width  ranges  from 
five  microns  to  forty  microns.  The  characteristics  of  the  transducers  and  the  phy¬ 
sical  aspects  of  the  pressure  distributions  will  be  discussed. 


Thursday ,  1  July  1993 
Conference  Room,  Bldg.  679 
Time:  10:30  AM 


POC:  Dr.  Promode  R.  Bandyopadhyay  (Code  8234;  x2588) 


1-IU/1-U2 
Reverse  Blank 


FLOWS  IN  MICRO  SYSTEMS 

Ho,  Chih-Ming 

CMS 

CAHl£L  fur  Mi£m  Systems 

University  of  California,  Los  Angeles 
Los  Angeles,  California  90024 
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OUTLINE 


•  Motivation 

•  Manufacturing  Process 

•  Experimental  Method 

•  Diffusion  and  Chemistry 

•  General  Etching  Model 

•  Pressure  Drop  in  Micro  Channel 
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1959 

THE  FIRST  PLANAR  TRANSISTOR-A 
TEARDROP  DEVICE  WHICH  MARKED  THE  BEGINNING  OF 
AN  INDUSTRY  REVOLUTION. 
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1964 

THE  FIRST  CONSUMER-ORIENTED  LINEAR  INTEGRATED 
CIRCUIT  IN  THE  WORLD. 
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1988 


MICRO  MOTOR 


We  are  witnessing  another  industrial  revolution. 


APPLICATIONS 


SENSOR  AND  ACTUATOR 


#  Micro  Prosthetic  Devices 

#  Biomedical  "Smart  Pills” 

#  Biological  Tissue  Connectors 

#  Biological  Micromanipulators 

#  Catheter-Based  Tools 

#  Microsurgical  Tools 

#  Clean  Room  Instrument 

#  Micro  Optical  Benches 

#  Micro  Light  Modulators 

#  Force-Balanced  Transducers 


#  Microvalves  and  Pumps 

#  Optical  Fiber  Connectors 

#  Microfilters 

#  Micro  Spacecraft 

$  Micro  Electrode  Arrays 
$  Electric  Switches  and  Relays 
$  Microelectromechanical  Muscle 

#  Electromechanical  Memories 

#  Miniature  Printer  Mechanisms 

#  Ultra-Precise  Positioning  Device 
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) 

SCIENCE  ISSUES 

•  Length  Scale 

•  Time  Scale 

•  Surface  Force  vs  Body  Force 

•  Material  Properties 
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□  □ 
□  □  □ 
D 


"Sacrificial"  Etching 


Sacrificial  layer: 
Structural  layer: 
Etchant: 
Selectivity: 


Silicon  Dioxide 
Polycrystaline  Silicon 
Hydrofluoric  Acid 
>  1000:1 
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POLYSILICON  MICROBEAM  TECHNOLOGY 


POLY -Si 


(R.T.  Howe  and  R.S.  Muller,  J.  Electrochem.  Soc.,  1983) 


Stator  Pole  Air  gap  Rotor  Bushing 


(b) 


(d) 


Fig.  6.  Micromotor  fabrication  process:  (a)  After  the  first 
LTO;  (b)  after  the  rotor-stator  polysilicon;  (c)  after  the  second 
LTO;  (d)  after  the  bearing  polysilicon;  and  (e)  released  motor. 
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Channels 

Width:  2jx  ~  200  p. 

Length:  3000  \l 

Thickness:  0. 2  jj,  ~  1.2  | x 

Acid 

HF:  1.8%  -49% 

14  concentrations 

Sacrificial  Material 

PSG 
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Experimental  Setup 


Microchannels 


Silicon  Nitride 


Microchannels 
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Governing  Equations 


0  8  (t)  x 

One-Dimentional  Eq.  of  Concentration 


dC(x,t)  3C(x,t)  02C(x,t) 

- 1-  u -  =  D - 

dt  dx  dx2 

B.C. 

c  [0,t]  =  Cb 

C[8(t),  t]  =  Cs  ??? 
Chemistry 

Si02  +  6HF  ^  H2SiF6  +  2H20 
Reaction  Flux 

Jr  (8)  =  f  [k,  CS(8)]  ??? 
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Convective  Term 


3C 

—  +  u 

at 


ac 

3x 


ac 

u  — 

3x 


u  ac/ax 

D  32c/3x2 
u  ??? 


U6 

D 
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Si02  +  6HF  -*  H2SiF6  +  2H20 

•  Assume  HF  &  H2SiF6 
dissolve  in  H20 

•  Change  of  the  volume  of  the 
solution  =  depletion  of  Si02 
+  production  of  H20 


•  Molecular  Weight  [MW] 


h2o 

18  gm/mole 

Si02 

60gm/mole 

•  Density 

h2o 

1  gm/cm3 

Si02 

2.1gm/cm3 

•  Depletion  of  1  mole  of 

Si02 

29cm3 

Gain  2  mole  of 

h2o 

36cm3 

•  Net  Volume  Gain  -  26% 
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J 


•  From  8  - 1  diagram 

— ►  d5/dt 

Vel.  of  solid/liquid  interface 

•  U  =  0.26d5  /dt 


u  dC/dx  ug 

D  32C/3x2  d 

~  [0.26  5  dS/dt]  /  D 
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Comparison  the  Convection  and  Diffusion  Terms  for  Different 

Concentration  of  HF 


8  (pm) 

■  49%  a  43.5%  ♦  38.2%  °  32%  1  26.5%  *  20%  •  13.8%  °  10%  x  7.1%  *  3.6%  +  1.8% 


Diffusion  Eq.: 


9t  3x2 

Diffusion  Time:  td  =  52(t)  /  D 


Etching  Time  :  t 

Ratio  between  etching  time  and 
diffusion  time: 

t*  =  t  /  td 
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Summary 


Small  convective  term 

ac  a2c 

- =  D  - 

at  ax2 

t  /  td  >>  1 


Concentraton 

=$  linear  distribution  in  x 


8(t)  jim 


Chemistry 


Fast  Reaction  -  Infinite  Reaction  Rate 


Cs  =  0 

Diffusion  Limited 
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Chemical  Reaction 


Reaction  Flux 


Jr(8)  =  f[k,  Cs(8)] 
Etching  Rate 

d8  /dt  ~  Jsi02 


dS  _  [MW]siq2 
dt  Psi02 


Jsi02 


_  3[MW]hf  Jhf 

PSi02  6 


~  f[k,  Cs(8)] 
Initial  Etching  Rate 

t  — )  0  8  — ) 


Cs(8)  — 3  Cb 

d8/dt]init.~  f[k,  Cb] 


Cs(5) 
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1000 


[Ui] 
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Chemical  Reactions 


HF  Equilibria : 


HF  v*  H+  +  F 


2  HF  ^H+  +  HF. 


iHFr±H+  +Hi-1F; 


Rate  Law :  ( Judge,  1971 ) 

-r  =  A[HF]  +  B[HF~]  +  C 

~  k1[HF]  +  k2[HF)2 
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Reaction  Flux  at  8 


Jr<8)  =  f[k,  Cs(8)]  =  kiCs(8)+k2Cs(8)2 
Diffusion  Flux  at  § 

Jd(8)  =  -D  dC(x)  /  dx 

=  D[Cb-Cs(5)]  /8  -  linear  distribution 

Jr(8)  =  Jd(5)  =  Jhf(5)  at  8  &  Eliminating  Cs 

Jhf(5)  =  2kJ  (  f  )2(1  +b(D)'<P) 

q>  =  Y  i  +  2b(^>  + 


where  b  =  kj  +  2Cbk2 
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d5  _  3[MW]hf  Jhf 

dt  pSi02  6 

=  a(  E-  )2(1  +  b(^)  -  q>) 
Solving  t  in  terms  of  5 


D{<p  <(>  +  2k^(|)+[kf(|)2-l]) 

8aCbkfk2(ki+Cbk2) 


D 

2ak? 


log[ 


<t>+ki<p 

2(k!+Cbk2) 


] 


where  <|>=  ki+2Cbk2+k?(-^-) 


a  =  {MW]hf  /  4psi02^2 


First-and-Second  Order  Model 


Assumptions : 

*  J[HF]  =  kl  Cs  +  k2  Cs2 

•  dC(x)  /  dx  =  [Cb-Cs(5>]  /8  (t) 

Results : 

D  =  1.6  x  10"5  cm^  /  s 
kj  =  1.2  x  10"4  cm  /  s 

V-2  -  6*5  x  10"^  cm4  /mole  s 

0  100  200 

Time  ( minute ) 
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Reaction  Flux  at  5 


Jr(8)  =  f[k,  Cs(8)]  =  kiCs(8)+k2Cs(8)2 
Diffusion  Flux  at  5 

Jd(§)  =  D[Cb-Cs(5)]  /8  -  linear  distribution 

Jr(8)  =  Jd(8)  at  8  &  Solving  Cs(S) 

Surface  Concentration 

Cg(S)  +  (  ^  )  Cs(8)  -  =  0 

^2  1^8  I128 
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Surface  Concentration  VS  Etch  Length 


PSG  Etching  Rate  ( cm/sec ) 


d8  /dt~  kiCs(S)+k2Cs(S)2 


SUMMARY 


•  Diffusion 

*  Convection  Effect  — *  0 

*  Concentration  Ihf~  x 

*  Asymptotic  d5/dt 

-  Diffusion  limited 

•  Chemistry 

*  kiCs(8)+k2Cs(5)2 

*  Initial  d8/dt 

-  Reaction  rate  limited 

•  Collapse  all  d8 /dt  with  Cs(8) 
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PRESSURE  DROP  IN  RECTANGULAR  CHANNELS 


•  Height:  h  =  1.2  pm 

•  Width:  w  =  40  Jim 

•  Length  =  2000  mm 

•  Distance  between  Pressure  Transducers:500  fun 

•  Piezoresistive  Transducer  :  250  jun  x  250  jun 

•  Transducer  Sensitivity  :  120n,V/(V*Psi) 
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Pressure  Sensors 


5 

MtCi 


Different!  Voltaic  (V) 
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Pressure  Distribution  Along  MicroChannel 


Pressure  Inside  Channel  (psi) 


PRVSVL5.XLS  Chart  21 


Comparison  of  Experimental  Results  and  Model  Prediction  for  Pressure  Drop 

Along  Channel 


Channel  length  (pm) 
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•  Kn  =  X  /  2h 


16p 

5pV2icRT 


•  Re  =  2h  U /  v 

•  M  =  U/  a 

•  Kn  ~  M  /  Re 


For  dp/dx  =  9  psi/pm 
V  =  28.  9  cm/sec 
Re  =  4.3  x  10  ‘2  M  = 

X  =  0.068  pm  Kn  : 


L-  =  8  V 

n  5V27I  hVRT 


V  =  10  V  • 
vHe  u  air 


7.7  x  10  -4 
=  0.028 
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Multi-disciplinary 

MANE,  EE,  Material  Science 


UCLA  &  Caltech 
8000  ft2  clean  room 

#  largest  one  in  U.S.  University 

#  compatible  for  micromachining 
&  microelectronics 


MICRO  SCALPEL 


Channel  Nozzle 


Micro  flow  meter 
Micro  pressure  transducer 
Microwave  wave  guide 
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Center  for  Micro.  Systems. 


CMS 


CMS  is  organized  for  the  purpose  of  advancing 
scientific  frontiers  and  manufacturing 
technologies  of  micro  systems.  This  center  will 
pursue  academic  excellence  in  fundamental 
research  and  will  actively  collaborate  with 
industries  and  government  for  technological 
developments. 
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Center  for  Micro.  Systems  CMS 

•  Science  =>  Technology  =>  Application 

•  15  Internationally  Renowned  Members 

•  Awarded  Research  Fund  -  US$  5.5  million 

•  Industrial  Participation 
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A  one-celled  protozoa,  a  euglena,  being  hold  by  the  microgripper.  The 
euglena  is  40  |im  long  and  7fim  in  diameter 
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JOURNAL  OF  MICROELECTr.OMECH.ASICAL  S'!  SI  FA'.S  VOL.  1.  NO.  1. 

1959 

There’s  Plenty  of  Room  at  the  Bottom 

Richard  P.  Feynman 


I  imagine  experimental  physicists  must  often  look  with 
envy  at  men  like  Kamerlingh  Onnes.  who  discovered  a 
held  like  low  temperature,  which  seems  to  be  bottomless 
and  in  which  one  can  go  down  and  down.  Such  a  man  is 
then  a  leader  and  has  some  temporary'  monopoly  in  a  sci¬ 
entific  adventure.  Percy  Bridgman,  in  designing  a  way  to 
obtain- higher  pressures,  opened  up  another  new  held  and 
was  able  to  move  into  it  and  to  lead  us  all  along.  The 
development  of  ever  higher  vacuum  was  a  continuing  de¬ 
velopment  of  the  same  kind. 

I  would  like  to  describe  a  held,  in  which  little  has  been 
done,  but  in  which  an  enormous  amount  can  be  done  in 
principle.  This  held  is  not  quite  the  same  as  the  others  in 
that  it  will  not  tell  us  much  of  fundamental  physics  (in  the 
sense  of,  “What  are  the  strange  panicles?”)  but  it  is  more 
like  solid-state  physics  in  the  sense  that  it  might  tell  us 
much  of  great  interest  about  the  strange  phenomena  that 
occur  in  complex  situations.  Furthermore,  a  point  that  is 
most  important  is  that  it  would  have  an  enormous  number 
of  technical  applications. 

What  I  want  to  talk  about  is  the  problem  of  manipulat¬ 
ing  and  controlling  things  on  a  small  scale. 

As  soon  as  I  mention  this,  people  tell  me  about  minia¬ 
turization,  and  how  far  it  has  progressed  today.  They  tell 
me  about  electric  motors  that  are  the  size  of  the  nail  on 
your  small  finger.  And  there  is  a  device  on  the  market, 
they  tell  me,  by  which  you  can -write  the  Lord's  Prayer 
on  the  head  of  a  pin.  But  that's  nothing;  that's  the  most 
primitive,  halting  step  in  the  direction  I  intend  to  discuss. 
It  is  a  staggeringly  small  world  that  is  below.  In  the  year 
2000,  when  they  look  back  at  this  age,  they  will  wonder 
why  it  was  not  until  the  year  1960  that  anybody  began 
seriously  to  move  in  this  direction. 

Why  cannot  we  write  the  entire  24  volumes  of  the  En¬ 
cyclopaedia  Britannica  on  the  head  of  a  pin  ? 

Let's  see  what  would  be  involved.  The  head  of  a  pin  is 
a  sixteenth  of  an  inch  across.  If  you  magnify  it  by  25  000 
diameters,  the  area  of  the  head  of  the  pin  is  then  equal  to 
the  area  of  all  the  pages  of  the  Encyclopaedia  Britannica. 
Therefore,  all  it  is  necessary  to  do  is  to  reduce  in  size  all 
the  writing  in  the  Encyclopaedia  by  25  000  times.  Is  that 
possible?  The  resolving  power  of  the  eye  is  about  1  /120 
of  an  inch— that  is  roughly  the  diameter  of  one  of  the  little 

MEMS  Editor's  Note:  This  manuscript  addresses  many  current  research 
issues.  It  is  the  transcript  of  a  talk  given  by  Richard  P.  Feynman  on  De¬ 
cember  26,  1939,  at  the  annual  meeting  of  the  American  Physical  Society 
at  the  California  Institute  of  Technology,  and  was  published  as  a  chapter 
in  the  Reinhold  Publishing  Corporation  book.  Miniaturization,  Horace  D. 
Gilbert,  Ed.  It  is  reprinted  with  the  consent  of  Van  Nostrand  Reinhold, 
New  York,  NY  10003. 

The  author,  deceased,  was  with  the  California  Institute  of  Technology. 
Pasadena,  CA. 
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dots  on  the  fine  half-tone  reproductions  in  the  Encyclo¬ 
paedia.  This,  when  you  demagnify  it  by  25  000  times,  is 
still  80  angstroms  in  diameter— 32  atoms  across,  in  an  or¬ 
dinary  metal.  In  other  words,  one  of  those  dots  still  would 
contain  in  its  area  1000  atoms.  So,  each  dot  can  easily  be 
adjusted  in  size  as  required  by  the  photoengraving,  and 
there  is  no  question  that  there  is  enough  room  on  the  head 
of  a  pin  to  put  all  of  the  Encyclopaedia  Britannica. 

Furthermore,  it  can  be  read  if  it  is  so  written.  Let's 
imagine  that  it  is  written  in  raised  letters  of  metal;  that  is, 
where  the  black  is  in  the  Encyclopaedia,  we  have  raised 
letters  of  metal  that  are  actually  1/25  000  of  their  ordi¬ 
nary  size.  How  would  we  read  it? 

If  we  had  something  written  in  such  a  way,  we  could 
read  it  using  techniques  in  common  use  today.  (They  will 
undoubtedly  find  a  better  way  when  we  do  actually  have 
it  written,  but  to  make  my  point  conservatively  I  shall  just 
take  techniques  we  know  today.)  We  would  press  the 
metal  into  a  plastic  material  and  make  a  mold  of  it,  then 
peel  the  plastic  off  very  carefully,  evaporate  silica  into  the 
plastic  to  get  a  very  thin  film,  then  shadow  it  by  evapo¬ 
rating  gold  at  an  angle  against  the  silica  so  that  all  the 
little  letters  will  appear  clearly,  dissolve  the  plastic  away 
from  the  silica  film,  and  then  look  through  it  with  an  elec¬ 
tron  microscope! 

There  is  no  question  that  if  the  thing  were  reduced  by 
25  000  times  in  the  form  of  raised  letters  on  the  pin,  it 
would  be  easy  for  us  to  read  it  today.  Furthermore,  there 
is  no  question  that  we  would  find  it  easy  to  make  copies 
of  the  master;  we  would  just  need  to  press  the  same  metal 
plate  again  into  plastic  and  we  would  have  another  copy. 


How  Do  We  Write  Small? 

The  next  question  is:  How  do  we  write  it?  We  have  no 
standard  technique  to  do  this  now.  But  let  me  argue  that 
it  is  not  as  difficult  as  it  first  appears  to  be.  We  can  reverse 
the  lenses  of  the  electron  microscope  in  order  to  demag¬ 
nify  as  well  as  magnify.  A  source  of  ions,  sent  through 
the  microscope  lenses  in  reverse,  could  be  focused  to  a 
very  small  spot.  We  could  write  with  that  spot  like  we 
write  in  a  TV  cathode  ray  oscilloscope,  by  going  across 
in  lines,  and  having  an  adjustment  which  determines  the 
amount  of  material  which  is  going  to  be  deposited  as  we 
scan  in  lines. 

This  method  might  be  very  slow  because  of  space 
charge  limitations.  There  will  be  more  rapid  methods.  We 
could  first  make,  perhaps  by  some  photo  process,  a  screen 
which  has  boles  in  it  in  the  form  of  the  letters.  Then  we 
would  strike  an  arc  behind  the  holes  and  draw  metallic 
ions  through  the  holes;  then  we  could  again  use  our  sys- 
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Recent  development  of  silicon  micromachining  technology  has  made  possi¬ 
ble  the  fabrication  of  micromechanical  parts  such  as  beams,  springs,  pin  joints, 
gears,  cranks,  sliders,  and  many  other  novel  structures  measured  in  microns. 
This  micromachining  technology  is  derived  and  modified  from  silicon 
integrated-circuit  technology.  Therefore,  it  also  has  many  important  advantages, 
such  as  precision  and  mass  production,  which  conventional  machining  tech¬ 
niques  can  not  offer.  Potential  applications  of  the  micromachining  technology 
are  many  and  wide.  In  fact,  this  new  field  of  micromachining  has  attracted 
attention  from  so  many  disciplines  including  physics,  chemistry,  engineering, 
biology,  medicine,  etc.  The  presentation  will  review  the  concept  and  technol¬ 
ogy  of  this  new  technology.  Examples  recently  developed  at  Caltech  will  be 
used  to  demonstrate  applications  of  micromachining  including  miniature  pres¬ 
sure  sensors,  Tera  Hz  wave  guides,  atomic  force  microscope  (AFM)  tips,  and 
neuron  prosthetic  devices.  Finally,  because  a  large  surface  to  volume  ratio  is  a 
distinctive  feature  of  these  microstructures,  many  microscience  issues  have 
become  extremely  important  in  the  fabrication,  characterization,  and  operation 
of  micromechanical  devices.  Among  the  issues,  micro  fluid  mechanics  and 
material  strength  will  specially  be  covered. 
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JOURNAL  OF  MICROELECTROMECHANICAL  SYSTEMS.  VOL.  I.  NO.  I,  MARCH  1992 


There’s  Plenty  of  Room  at  the  Bottom 

Richard  P.  Feynman 


I  imagine  experimental  physicists  must  often  look  with 
envy  at  men  like  Kamerlingh  Onnes,  who  discovered  a 
field  like  low  temperature,  which  seems  to  be  bottomless 
and  in  which  one  can  go  down  and  down.  Such  a  man  is 
then  a  leader  and  has  some  temporary  monopoly  in  a  sci¬ 
entific  adventure.  Percy  Bridgman,  in  designing  a  way  to 
obtain' higher  pressures,  opened  up  another  new  field  and 
was  able  to  move  into  it  and  to  lead  us  all  along.  The 
development  of  ever  higher  vacuum  was  a  continuing  de¬ 
velopment  of  the  same  kind. 

1  would  like  to  describe  a  field,  in  which  little  has  been 
done,  but  in  which  an  enormous  amount  can  be  done  in 
principle.  This  field  is  not  quite  the  same  as  the  others  in 
that  it  will  not  tell  us  much  of  fundamental  physics  (in  the 
sense  of,  “What  are  the  strange  particles?1  *)  but  it  is  more 
like  solid-state  physics  in  the  sense  that  it  might  tell  us 
much  of  great  interest  about  the  strange  phenomena  that 
occur  in  complex  situations.  Furthermore,  a  point  that  is 
most  important  is  that  it  would  have  an  enormous  number 
of  technical  applications. 

What  I  want  to  talk  about  is  the  problem  of  manipulat¬ 
ing  and  controlling  things  on  a  small  scale. 

As  soon  as  I  mention  this,  people  tell  me  about  minia¬ 
turization,  and  how  far  it  has  progressed  today.  They  tell 
me  about  electric  motors  that  are  the  size  of  the  nail  on 
your  small  finger.  And  there  is  a  device  on  the  market, 
they  tell  me,  by  which  you  can  write  the  Lord's  Prayer 
on  the  head  of  a  pin.  But  that's  nothing;  that's  the  most 
primitive,  halting  step  in  the  direction  I  intend  to  discuss. 
It  is  a  staggeringly  small  world  that  is  below.  In  the  year 
2000,  when  they  look  back  at  this  age,  they  will  wonder 
why  it  was  not  until  the  year  1960  that  anybody  began 
seriously  to  move  in  this  direction. 

Why  cannot  we  write  the  entire  24  volumes  of  the  En¬ 
cyclopaedia  Britannica  on  the  head  of  a  pin  ? 

Let's  see  what  would  be  involved.  The  head  of  a  pin  is 
a  sixteenth  of  an  inch  across.  If  you  magnify  it  by  25  000 
diameters,  the  area  of  the  head  of  the  pin  is  then  equal  to 
the  area  of  all  the  pages  of  the  Encyclopaedia  Britannica. 
Therefore,  all  it  is  necessary  to  do  is  to  reduce  in  size  all 
the  writing  in  the  Encyclopaedia  by  25  000  times.  Is  that 
possible?  The  resolving  power  of  the  eye  is  about  1  / 120 
of  an  inch — that  is  roughly  the  diameter  of  one  of  the  little 

AtEMS  Editor's  Note:  This  manuscript  addresses  many  current  research 
issues.  It  is  the  transcript  of  a  talk  given  by  Richard  P.  Feynman  on  De¬ 
cember  26,  1959,  at  the  annual  meeting  of  the  American  Physical  Society 
at  the  California  Institute  of  Technology,  and  was  published  as  a  chapter 
in  the  Reinhold  Publishing  Corporation  book.  Miniaturization,  Horace  D. 
Gilbert,  Ed.  It  is  reprinted  with  the  consent  of  Van  Nostrand  Reinhold, 
New  York,  NY  10003. 

The  author,  deceased,  was  with  the  California  Institute  of  Technology, 
Pasadena,  CA. 
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dots  on  the  fine  half-tone  reproductions  in  the  Encyclo¬ 
paedia.  This,  when  you  demagnify  it  by  25  000  times,  is 
still  80  angstroms  in  diameter— 32  atoms  across,  in  an  or- 
dinaiy  metal.  In  other  words,  one  of  those  dots  still  would 
contain  in  its  area  1000  atoms.  So,  each  dot  can  easily  be 
adjusted  in  size  as  required  by  the  photoengraving,  and 
there  is  no  question  that  there  is  enough  room  on  the  head 
of  a  pin  to  put  all  of  the  Encyclopaedia  Britannica. 

Furthermore,  it  can  be  read  if  it  is  so  written.  Let's 
imagine  that  it  is  written  in  raised  letters  of  metal;  that  is, 
where  the  black  is  in  the  Encyclopaedia,  we  have  raised 
letters  of  metal  that  are  actually  1  /2 5  000  of  their  ordi¬ 
nary  size.  How  would  we  read  it? 

If  we  had  something  written  in  such  a  way,  we  could 
read  it  using  techniques  in  common  use  today.  (They  will 
undoubtedly  find  a  better  way  when  we  do  actually  have 
it  written,  but  to  make  my  point  conservatively  I  shall  just 
take  techniques  we  know  today.)  We  would  press  the 
metal  into  a  plastic  material  and  make  a  mold  of  it,  then 
peel  the  plastic  off  veiy  carefully,  evaporate  silica  into  the 
plastic  to  get  a  very  thin  film,  then  shadow  it  by  evapo¬ 
rating  gold  at  an  angle  against  the  silica  so  that  all  the 
little  letters  will  appear  clearly,  dissolve  the  plastic  away 
from  the  silica  film,  and  then  look  through  it  with  an  elec¬ 
tron  microscope! 

There  is  no  question  that  if  the  thing  were  reduced  by 
25  000  times  in  the  form  of  raised  letters  on  the  pin,  it 
would  be  easy  for  us  to  read  it  today.  Furthermore,  there 
is  no  question  that  we  would  find  it  easy  to  make  copies 
of  the  master,  we  would  just  need  to  press  the  same  metal 
plate  again  into  plastic  and  we  would  have  another  copy. 


How  Do  We  Write  Small? 

The  next  question  is:  How  do  we  write  it?  We  have  no 
standard  technique  to  do  this  now.  But  let  me  argue  that 
it  is  not  as  difficult  as  it  first  appears  to  be.  We  can  reverse 
the  lenses  of  the  electron  microscope  in  order  to  demag¬ 
nify  as  well  as  magnify.  A  source  of  ions,  sent  through 
the  microscope  lenses  in  reverse,  could  be  focused  to  a 
very  small  spot.  We  could  write  with  that  spot  like  we 
write  in  a  TV  cathode  ray  oscilloscope,  by  going  across 
in  lines,  and  having  an  adjustment  which  determines  the 
amount  of  material  which  is  going  to  be  deposited  as  we 
scan  in  lines. 

This  method  might  be  very  slow  because  of  space 
charge  limitations.  There  will  be  more  rapid  methods.  We 
could  first  make,  perhaps  by  some  photo  process,  a  screen 
which  has  holes  in  it  in  the  form  of  the  letters.  Then  we 
would  strike  an  arc  behind  the  holes  and  draw  metallic 
ions  through  the  holes;  then  we  could  again  use  our  sys- 
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SILICON  MICROMACHINING 

•  BULK  MICROMACHINING  -  SUBSTRATE 

—  Isotropic  and  anisotropic  etchants  (KOH,  NaOH,  EDP) 

—  Etching  stop  (p+  boron) 


•  SURFACE  MICROMACHINING  -  THIN  FILMS 

—  Structural  materials  (Si,  S^N^  metals) 

—  Sacrificial  layers  (Si02,  Photoresist) 

—  Selective  etchants 


<1 10>  Surface  Orientation 
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COMPARISON  OF  ANISOTROPIC  ETCHANTS  FOR  SILICON  [5.8,5.9,5.10] 


ETCHANT 

TYPICAL 

COMP. 

TEMP. 

(Diluent) 

°C 

KOH 

44  gr 

85 

(Water) 

100  ml 

NaOH 

10  gr 

65 

(Water) 

100  m! 

Ethylene 

diamine 

750  ml 

Pyrocatecho! 

120  gr 

115 

(Water) 

240  ml 

h2n4 

100  ml 

100 

(Water) 

100  ml 

(100) 

ETCH 

RATE 

p.m/mln 

(100)/(111) 
ETCH  RATE 
RATIO 

MASK 

(ETCH  RATE) 

1.4 

400:1 

SI3N4 

SIO2  (14  A/mln) 

0.25-1.0 

>300:1 

SI3N4 

SIO2  (7  A/mln) 

1.25 

35:1 

SIO2  (2  A/mln) 
SI3N4  (1  A/mln) 
Ta'AU'Cr.Ag'Cu 

2.0 

10:1 

SIO2  (<2  A/mln) 

0  I.  SMrru 


well  500x500  pn 
depth  180  pn 


S.E.I 


na9nified  tilted 
vieu 


00001  lOOpn 


f^^cr op es  5^^00pn| 

no  (>  n  i  no  1  0  un 

*  \  •'I 

k/  tilkZs.E.k/ 

h/  U  \ 

kZ  I 

kZ  LZ  kZ-j 

1/  ZZ  ! 

!ZZ.  1 

kz  kZ  kZ 

[v  ! 

h/. 

kZ  kZ  kz 

\/ 

!kZ 

kZ  .kZ  kZ1 

t.'r- 


00009  lOOOpn 
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HORIZONTAL  UMOST  VERTICAL  UMOST 

Figure  5.  Two  MOS  type  devices  on  the  (110).  The  dotted  lines 
show  the  inversion  layers. 

D.L.  Kendall  &  G.R.  de  Guel,  Workshop  on  M.  M.  T.,  Cleveland,  1984. 


SN 


Figure  6. 

a.  Intercalated  studs  for  high  area 
cathode  (13). 

b.  Tall  studs  showing  the  effect  of 
a  dislocation. 

c.  MJ-STRATE  showing  vents,  conduc¬ 
ting  bus-lines,  logic  intercon¬ 
nects  on  back  side  and  active  de¬ 
vices  on  top. 

d.  Well  aligned  grooves  on  the  (231) 


,+ 


CONCENTRATION/ORIENTATION 
DEPENDANT  ETCHANT 

Example:  KOH,  Isopropyl,  H20-  R  >60R 

K  *  (IOO)  (III) 

EDP-R(|00)>50R(|||);  Etch  stops  if 
Na  I0,9crrr3 
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Fig,  61  Relative  silicon  etch  rate  as 
a  function  of  boron  dopant  concentra¬ 
tion  for  EDP  type  S, 
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POLYSILICON  MICROBEAM  TECHNOLOGY 


POLY -Si 


(R.T.  Howe  and  R.S.  Muller,  J.  Electrochem.  Soc 1983) 


SURFACE  MICROMACHINING 


#  Structural  Materials 

#  Sacrificial  Materials 

#  Selective  Etchants 
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SEM  photographs  of  the  cantilever  and  spring  joint  shutters 


Spring  Joint 
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mu  lJb  lv  Mb  _  ___  - 

Mfc  ttUBHffW  Tut  Q:  1ST  S*ll  H;  WL  iWB 


8S  litres 
tail:  8, 

Mil:  211 .525  Cicrest 


tall:  SJBv  !iij: 


'  rrn\ 
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ICil  II:  IT 
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Mil:  9.56  Kv  Mag:  0.78  Kx 

lift:  9.584  M I  crons  Test  ID:  U-MECH  Sanple  ID:  SPRING 
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Caltech  Micromachining  Lab 

Space 

#  Surface  micromachining  area  (class  100,  800  sq.  ft.) 

#  Bulk  micromachining  area  (class  10,000,  700  sq.  ft.) 

#  CAD  and  measurement  area  (500  sq.  ft.) 

Equipment 

#  GCA  stepper  and  contact  aligner 

#  LPCVD  polysilicon  reactor 

#  LPCVD  Si3N4  (SiN)  reactor 

#  LPCVD  Si02  (PSG)  reactor 

#  Diffusion  furnaces 

#  PECVD  reactor 

#  Plasma  etcher 

#  RIE  etcher 

#  Evaporators 

Planned  Equipment 

#  LPCVD  Tungsten 

#  Electroforming  Tank 

#  Sputtering  machine 
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CURRENT  CALTECH  PROJECTS 


#  Culture  Neuron  Probes 

#  Si  STM  and  AFM  Probes 

#  Integrated  Tera  Hz  Waveguide  Circuits 

#  Miniature  Disc  Head  (Micromagnetics) 

#  Hot-point  Anemometer 

#  Surface-Micromachined  Pressure  Sensors 

#  MicroChannel  Flow 

#  Accelerometer 
.  #  Micro  Gyros 

#  Micromechanical  Properties 
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Figure  1  A  culture  dish  mounted  on  a  Figure  2  The  central  electrode  pattern  of 

glass  substrate.  the  dish  in  Fig.  1. 


Insulated  gold  wire 
to  electronics 


r 


1-196 


1-197 


NEUROCHIP  FABRICATION  PROCESS 


□  -Si 


-  metal 


-  insulation 


□  -  Si02 


j  -  B  doped  Si 


-  Si3N4 
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Itt  Vitro  Use  of  Neurochips 


embryonic  cell 


a)  Neuron  after  implantation 


fju^ro^  ^neurite 


gold  lead  platinum  black 

b)  Neuron  after  growth 
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A  Live  Neural  Network 
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} 

NEUROPROBE 


SILICON  NEUROPROBE 

NEURON  WELLS,  NO  ELECTRODES 


SILICON  NEUROCHIP 

NEURON  WELLS  WITH  ELECTRODES 


FOR  IN  VIVO  STUDIES,  CENTRAL  FOR  IN  VITRO  STUDIES;  NEURON 

NERVOUS  SYSTEM;  NEURON  SURVIVABILITY  AND  EXCITABILITY 

SURVIVABILITY  AND  GROWTH 


SILICON  NEUROPROBE 

NEURON  WELLS  WITH  ELECTRODES 

FOR  IN  VIVO  STUDIES; 
STIMULATION  AND  RECORDING 
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1  ASI  K 


AFM  Imaging 


Non-Contact  Modes 


Contact  Modes 
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IODiD  u  t— H 
DS-I  00>D  OS  005  010 
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55-100:0  00  005  Dll 
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C  altech  Micromachining  Laboratory 


ret  Propulsion  Laooraton. 


WR-10  WAVEGUIDE  FABRICATION  PROCESS 


KOH  etching 


(110)  silicon  (50  mils) 


silicon  mtnde 


(110)  silicon  (50  mils; 


boron  oxide/polyimide 
bonding 


(100)  silicon  (50  mils) 


chrome  (250  A 
gold  (5000  A* 
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IMPORTANT  MILLIMETER  WAVE  AND  SUBMILLIMETER  WAVE  APPLICATIONS 


0  Mixer  and  Frequency  Multiplier  Embedding  Circuits 

0  Directional  Couplers 

0  Waveguide  Transformers 

°  Waveguide  -  to  -  Planar  Circuit  Transitions 

°  Low  Loss  Filters 

0  Rectangular  and  Conical  Feedhorns 

0  Dichroic  Filter  Plates 


Si  Micromachlned  waveguide  C  frcuit 


Ijx  SiN  Membrane 


Active  Device 


Si'icor.  wafer  with  several 
etched  waveguides  and 
devices 


Etched 

Waveguide 


RF  Filter 


RF  in 


IF,  LC 
Out 
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Current  Research  Activities 
at  UCLA/CalTech 


(1)  supported  by  AMC,  HP,  IBM,  Maxtor,  Quantum,  Read-Rite  and  Seagate  Tfechnology  Corporations 
with  matching  funds  from  the  State  of  California  MICRO  program. 

^  supported  by  the  National  Storage  Industry  Consortium  (NSIC). 

(3)  proposal  submitted  to  the  Defense  Advanced  Research  Program  Agency  (DARPA). 
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Form  Factor  Food-Chain 


Palmtop 


Notebook 


Laptop 


SCALED  Technology  -  Super  Compact  Array  of 
Low-cost  Enormous-capacity  Disk-drives 


Integrated  Read-Write-Head/Suspension/ 
T  orsional-Bearings/V  oice-Coil-Motor 

(Long-Term  Goal) 


Silicon  Micromachined  Gimbal/Suspension 

with  integrated  Electrodes 


* 


100  iO  u 

ntJ-l  M,0  PD  m  OOP 
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Some  Recently  Fabricated  Silicon 
Microgimbals  (beam  width  ~  80  microns) 


ilH.MV  '♦Mill, 


069305  25KV  X36.0  830um 
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HOT-POINT  ANEMOMETER 


Aluminum 


diffusion 

barrier 


sensing 

resistor 


Si  substrate 
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V-l  Characteristic 
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■  Caltech  Micromachining  Lab  ■ 


Flow  Channel  &  Pressure  Sensors  Design 


Gas  Inlet 


500  urn 


Bridging 

Channel 

/  Outlet 


Pressure  Sensors 


( 250x250  fim ) 
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■  Caltech  Micromachining  Lab  ■ 


Pressure  Sensor 

Piezoresistors 
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■  Caltech 


Etch  Window 


Lab  ■ 


lire  Sensor 


Flow  Channel 


Bridging  Channel 
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■  Caltech  Micromachining  Lab  ■ 


Micro  Flow  System  Chip 

Flow  Channels 


(250x250nm2) 
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■  Caltech  Micromachining  Lab  ■ 


Fabrication  Process 
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Pressure  Sensor  for  Integrated  Microchannel/Pressure  Sensor  System 


SEM  photo  of  a  pressure  sensor  Pressure  sensor  calibration  curve 

(Pressure  applied  from  outside  of 
the  sensor.) 


Differential  Voltage  ( V ) 


■  Caltech  Micromachining  Lab  ■ 


Experimental  Results  (  Pressure  Sensor ) 


0  5  10  15  20 

Pressure  ( psig ) 
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■  Caltech  Micromachining  Lab  ■ 


Experimental  Results  (  Micro  Channel ) 


p(x)/p(0) 


x  (pm) 


Pressure  distribution  of  air  flow  in  a  microchannel  under  various 
pressure,  D  =  1  um,  L  =  2000  urn. 


) 


Figure  5 
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V(x)  (m/sec) 


ORIFICE  TUBE 


VISCOUS  DRAG  TUBE 


Fig.  4:  Mean  Flow  Rate  Sensors 
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Flow/Sensor  Chip  ( Next  Generation ) 


Inlet/Outlet 


Thermistors 


HP 

uu 
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Pressure  Sensors 


2D  Micro  Jets 
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POTENTIAL  MICROFLOW  APPLICATIONS 


1.  Microsurgical  tools 

2.  Micromanipulators 

3.  Air  bearing 

4.  Microfluidics 


5.  Chemical  sensors 

6.  Flow  sensors 

7.  Chromatograph 

8.  Biomedical  devices 
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Typical  Applications  for  Micromechanical  Structures  in  Silicon 


Cryogenic  Microconnectors 

Microprobes 

Fiber  Optic  Couplers 

Micropumps 

Film  Stress  Measurement 

Microswitches 

Fluidic  Components 

Microvacuum  Tubes 

1C  Coolers 

Microvalves 

Ink  Jet  Nozzles 

Nerve  Regenerators 

Laser  Beam  Deflectors 

Photolithography  Masks 

Laser  Resonators 

Pressure  Switches 

Light  Modulators 

Pressure  Regulators 

Membranes 

Programmable 

Microaligners 

RMS  Converters 

Microbalances 

Thermal  Print  Heads 

Microfuses 

Thermopile 

Microgears 

Torsional  Mirrors 

Micromolds 

Vibrating  Microstructures 

Micromotors 

Wind  Tunnel  Skin  Friction  Gages 

Micropositioners 
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SENSOR  APPLICATIONS 

1.  Flow  Sensors  5.  Force  Sensor 

2.  Humidity  Sensors  6.  Pressure  Sensors 

3.  Chemical  Sensors  7.  Acceleration  Sensors 

4.  Infrared  Sensors  8.  Tactile  Sensors 
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LONG  TERM  INTERESTS 


***  1.  Microelectronics  + 

\ 

\ 

'Microsensors  + 

.5 

Micromachines  (sM's). 

2.  Implantable  “Insulin  Pill.” 

3.  Micro-gyroscope. 

4.  Micro-robots  and  Micro-spacecraft. 
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Small  Machines, 

,  Large  Opportunities: 

i  A  Report  on  the  Emerging 
j  Field  of  Microdynamics 

i - 

Report  of  the  NSF  Workshop 
on  Microelectromechanical 
!  Systems  Research 


POSSIBLE  APPLICATIONS 


#  Micro  Prosthetic  Devices 

#  Biomedical  “Smart  Pills" 

#  Biological  Tissue  Connectors 

#  Biological  Micromanipulators 

#  Catheter-Based  Tools 

#  Microsurgical  Tools 

#  Clean  Room  Instrument 

#  Micro  Optical  Benches 

#  Micro  Light  Modulators 

#  Force- Bala  need  Transdu 


#  Microvalves  and  Pumps 

#  Optical  Fiber  Connectors 

#  Microfilters 

#  Micro  Spacecraft 

#  Micro  Electrode  Arrays 

#  Electric  Switches  and  Relays 

#  Microelectromechanical  Muscle 

#  Electromechanical  Memories 

#  Miniature  Printer  Mechanisms 

#  Ultra-Precise  Positioning  Device 
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In-Situ  Monitoring  and  Universal  Modelling  of  sacrificial 
SiC>2  Etching  by  Hydrofluoric  Acid 

Jianqiang  Liu,  Yu-Chong  Tai 

MS  116-81,  Caltech,  Pasadena,  CA  91125 

Kin-Cheok  Pong,  Jiajing  Lee,  Yitshak  Zohar,  Chin-Ming  Hu 

MANE,  UCLA,  Los  Angelas,  CA  90024 


Summary 

A  video  system  is  designed  to  do  in-sit u  and  accurate 
monitoring  of  sacrificial  phosphosilicate  glass  (PSG)  etch¬ 
ing  in  microchannels  by  hydrofluoric  acid.  New  phenomena 
including  bubble  formation  and  liquid  condensation  have 
been  found  to  either  enhance  or  hinder  etching  process  de¬ 
pending  on  the  size  and  geometry  of  the  channels.  For 
the  first  time,  an  universal  model  for  microchannel  etch¬ 
ing  is  identified  to  cover  a  wide  range  (10-49  wt%)  of  HF 
concentration.  This  model  includes  chemical  diffusion  and 
second-order  chemical  reaction. 

Abstract 

Surface  micromachining  requires  chemical  etching  of  sac¬ 
rificial  layers.  Silicon  dioxide  or  PSG  etching  in  HF-based 
solution  is  one  of  the  most  frequently  used  techniques  [1,2]. 
In  most  of  the  surface  micromachining,  the  length  of  the 
'erificial  layer  is  short  (  <  100  /*m)  so  that  there  are  no 
^ous  problems  such  as  unbearable  long  etching  time  and 
structural  damage.  However,  several  Caltech  projects  in¬ 
volve  large  silicon-nitride  diaphragms  and  long  microchan¬ 
nels,  and  prolonged  exposure  of  these  structures  to  HF  (> 
40  wt%)  greatly  reduces  yield.  To  solve  the  problems,  one 
needs  a  mode]  that  covers  a  wide  range  of  HF  concentra¬ 
tion,  hence  called  a  universal  model.  This  is  because  the 
starting  HF  concentration  could  be  different,  and,  once 
etching  starts,  chemical  concentration  inside  the  sacrificial 
gap  during  etching  varies  as  a  function  of  time  and  position. 
Many  works  have  been  done  on  modelling  silicon-dioxide 
etching,  e.g.  Monk  et.  al.[3]  and  Helms  et.  al.  [4];  How¬ 
ever,  no  universal  model  has  been  identified  so  far.  In  the 
following,  we  describe  our  effort  in  search  of  such  a  model. 

First,  a  video  system  is  designed  to  do  in-situ  and  accu¬ 
rate  monitoring  of  sacrificial  phosphosilicate  glass  (PSG) 
etching  (Fig.  1).  This  system  basically  consists  of  a  micro¬ 
scope,  a  video  camera  and  a  tape  recorder.  A  polypropy¬ 
lene  container  covered  by  a  transparent  plastic  slide  is  used 
to  take  sample  dies  and  HF.  For  etching  at  different  tem¬ 
peratures,  a  hot  chuck  is  used  underneath  the  container 
with  calibrated  temperature  control.  The  container  can 
take  a  1  by  1  cm2  silicon  die,  which  has  various  straight 
PSG  (~  6wt  %)  channels  covered  by  LPCVD  low-stress 
'icon  nitride  (Fig.  2).  Each  channel  has  a  ruler  nearby 
Ah  10  micron  divisions  for  convenient  reading  of  the  etch¬ 
ing  length.  Many  experiments  have  been  done,  and  a  large 
amount  of  repeatable  data  (  it  5  %)  have  been  collected  at 
different  HF  concentrations  and  temperatures  as  in  Figs. 
3  and  4. 

New  phenomena  have  been  found  during  the  etching  ob¬ 
servation.  For  example,  bubbles  inside  the  channels  may 


form  during  the  HF  etching.  Interestingly,  when  channel 
width  is  larger  than  a  threshold  value  (  ~  4  /im),  bubbles 
are  movable.  They  always  move  from  the  nucleation  site  to 
the  etching  window.  In  this  case,  bubbles  tend  to  enhance 
etching  by  probably  introducing  micro-flows  in  the  channel. 
On  the  other  hand,  when  width  is  below  the  threshold  value 
bubbles,  if  formed,  become  stationary  but  grow  big  enough 
to  block  the  channel.  These  bubbles  then  significantly  hin¬ 
der  the  etching  of  the  channels.  In  general,  the  formation 
of  bubble  doesn’t  always  happen,  and  it  depends  on  the 
temperature  and  the  quality  of  PSG.  It  is  also  found  that 
after  the  bubble  formation  liquid  condensation  can  happen 
on  the  silicon  nitride  surface  containing  the  bubbles. 

Our  data  are  collected  from  video  tape  with  an  accurate 
enough  time  resolution  of  33  ms  to  examine  many  possible 
models.  The  result  is  summarized  in  Table  1.  Pure  diffu¬ 
sion  model  doesn’t  fit  the  data  well  enough.  Deal-Grove 
model  (diffusion  and  first  order  chemical  reaction)  can  only 
fit  data  with  a  fixed  HF  concentration.  In  other  words,  it 
doesn’t  cover  a  wide  range  of  HF  concentration.  Langmuir- 
Hinshelwood  model  (spedal  chemical  kinetics  coupled  with 
diffusion)  is  similar  to  Deal-Grove  model.  It  turns  out 
that  only  a  model  using  diffusion  and  second-order  chem¬ 
ical  reaction  (r=-kC2WFj)  universally  fit  all  data  covering 
a  wide  HF  concentration  range  from  10-49  wt  %  with  a 
k=0.1  cml/molesec.  and  diffusion  constant  D=1.5  x  10~6 
cm2/sec  (Fig.  3).  The  accuracy  of  this  model  is  about 
10%. 

Finally,  temperature  effects  on  sacrificial  layer  etching 
have  been  also  studied(Fig.  4).  It  is  found  that  chemical 
reaction  constant  is  affected  more  than  that  in  diffusion 
constant.  Including  the  temperature  effects  in  the  model 
has  also  been  done.  A  lot  more  of  the  details  will  be  pre¬ 
sented  if  this  abstract  is  accepted. 
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Etched  Channel  l^enfth  (micron) 


PSG  Etching  window 


Silicon  nitride 


Figure  1.  Schematic  of  experimental  setup.  Vapor  condensation 
is  avoided  by  using  a  plastic  cover  on  the  HF  container. 
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Figure  3.  Typical  PSG  channel  etching  data.  Solid  lines  are  modelled 
using  diffusion  and  a  second  order  chemical  reaction  and  diffusion 
model  fit. 


Silicon  sbustrate 


Figure  2.  (a)  An  etched  microchannel  with  a  ruler.  The  silicon 
nitride  diaphgram  was  flat  during  etching,  (b)  Schematic  of  its 
cross-section. 
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Firgure  4.  Temperature  and  bubble  effects  on  channel  etching. 


Table  1 .  Model  List 


Model 


Parameters 


Fitting 


Chemical 
Rate  Law 


Diffusion 


Inadequate 


infinite 
reaction  rate 


Deal -Grove 


Inadequate 


r=-k[HF] 


D,  k 


Langmuir- 

Hinshelwoo 


r=  -kl  [HF]/ 
(l+k2[HF]) 


D,  kljc2  Inadequate 


Second-order  i^-k[HF]*n  D,  k,  n 


Good,  n=2 


Figure  5.  Bubble  formation  and  its  sequencial  movement,  (a)  A 
bubble  forms  at  the  interface  between  HF  and  PSG.  (b)  The  bubble 
expands  and  moves  toward  etching  window,  (c)  The  bubble 
saturates  at  etching  window.  Black  circular  spots  are  accumlaied 
eas  around  windows. 
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A  MICROMACHINED  SILICON  STRUCTURE 
FOR  IN  VITRO  STUDIES  OF  NEURAL  NETWORKS 

Svetlana  Tatic-Lucic*,  Yu-Chong  Tai,  John  Wright 
Department  of  Electrical  Engineering,  California  Institute  of  Technology,  Pasadena,  CA  91125,  U.S.A 

Jerome  Pine,  Timothy  Denison 

Department  of  Biology,  California  Institute  of  Technology,  Pasadena,  CA  91125,  U.S.A 


Summary — This  paper  describes  a  novel  structure  for  an  extracellular  recording.  A  4x4  array  of  neuron 
wells  has  been  fabricated  on  a  20  /im  thick  silicon  membrane  using  double  side  photolitography.  The  neuron 
well  is  25  pm  square,  with  a  circular  gold  electrode  3  pm  diameter  at  the  bottom  and  the  heavily  doped 
silicon  grillwork  at  the  top.  The  experiments  confirmed  that  this  structure  is  biocompatible.  In  this  way 
simultaneous  recording  with  one-to-one  correspondence  between  the  electrode  and  the  cell  can  be  performed 
from  the  sixteen  neurons. 

Abstract — The  understanding  of  the  functioning  of  the  nervous  system  became  an  imperative  in  the 
scientific  world.  One  way  to  achieve  that  difficult  task  is  to  do  studies  on  so-called  cultured  neurons.  They  are 
the  individual  neurons,  isolated  from  the  brain  and  grown  in  the  solution  which  has  all  neccessary  components 
essential  for  their  growth  (sugars,  proteins, etc.)  The  goal  is  to  measure  the  electric  signals  originating 
from  the  cultured  neurons.  Different  techniques  were  developed,  both  intracellular  and  extracellular.  The 
extracellular  techniques  are  not  invasive  or  harmful  for  the  cell,  although  the  signal-to-noise  ratio  is  not 
as  good  as  with  intracellular  recording.  Easy  fabrication  of  the  silicon  microdevices  has  a  potential  for  the 
significant  improvement  of  the  existing  methods. 

The  ancestors  of  our  structure  (in  further  text,  we’ll  refer  to  it  as  a  ’’neurochip”)  are  the  diving  board 
electrode  (Figl(a))  and  the  polyimide  wells  (Figl(b)),  both  previously  fabricated  at  CALTECH.  They  en¬ 
abled  the  long  term  recording  and  stimulation  of  the  neurons,  but  the  diving  board  electrodes  were  difficult 
to  mount  and  unable  to  detect  postsynaptic  potentials,  and  the  polyimide  wells  were  not  suitable  for  the 
mass  production. 

The  neurochip  represents  the  combination  of  the  well  structure  with  an  electrode  structure  like  that  of 
the  diving  board.  It  has  a  4x4  array  of  the  silicon  wells  situated  at  the  center  of  a  9x3  mm  silicon  membrane, 
20  /xm  thick.  The  structure  of  the  individual  neuron  well  as  well  as  the  dimensions  and  the  design  of  the 
neurochip  are  shown  in  Fig  2. 

The  fabrication  process  begins  with  formation  of  6  pm  diameter  nitride  dots,  which  determine  the 
position  of  the  neuron  wells  electrodes.  The  LOCOS  process  is  the  used  to  create  a  1  pm  oxide  step  in 
the  bottom  of  the  neuron  well.  Then  a  Cr/Au/Cr  metalization  is  deposited  using  a  lift-off  process.  The 
insulation  layer  is  deposited  after  the  metalization.  The  holes  in  dielectric  are  opened  to  expose  the  bonding 
pads.  The  windows  are  etched  in  the  insulation  layer  at  the  back  and  succeeding  anisotropic  EDP  etch 
removes  silicon  all  the  way  to  the  heavily  boron-doped  layer.  Then  the  grillwork  for  the  neuron  wells  is 
patterned  at  the  bottom  of  the  cavity,  using  the  projection  stepper.  R1E  etching  selectively  removes  the 
heavily  doped  boron  layer  and  the  neuron  wells  are  then  formed  with  am  EDP  etch. 

Fig.3  shows  the  SEM  of  the  single  neuron  well.  The  small  circle  at  the  bottom  of  the  well  is  a  gold 
electrode.  The  network  of  the  successfully  growing  rat  superior  cervical  ganglion  (SCG)  neurons  is  shown  at 
Fig. 4.  The  picture  shows  the  neurites  growing  out  of  the  wells  from  the  invisible  neurons  implanted  in  them 
and  connecting  among  themselves. 

*  The  author  to  whom  correspondence  should  be  addressed. 

Phone:  (818)356-3884 

Fax:  (818)356-2944 
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Figure  2:  a)  Neuron  just  implanted  in  a  well,  b)  Neuron 
increases  in  size  and  grows  processes  out  of  the  well. 
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Figure  4:  SEM  of  a  single  neuron  well  25x25/on2  at  the  top 
with  3//m  circular  gold  electrode  at  the  bottom. 


Figure  5:  A  live  neural  network  formed  by  rat  superior  cer¬ 
vical  ganglion  neurons.  Clearly  seen  are  the  processes  grow¬ 
ing  out  of  the  wells. 

1-254 


Figure  3:  SEM  of  a  4x4  array  of  neuron  wells  at  the  bottom 
of  a  neurochip. 
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Micromachined  Channel/Pressure  Sensor  Systems  for  Micro 

Flow  Studies 

Jianqiang  Liu  &  Yu-Chong  Tai 
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Kim-Cheok  Pong  &  Chih-Ming  Ho 

MANE,  UCLA,  Los  Angeles,  CA  90021 


Summary 

For  the  first  time,  integrated  microchannel /pressure  sensor 
systems  have  been  designed  and  fabricated  for  direct  mea¬ 
surement  of  pressure  distribution  along  a  microchannel.  Both 
bulk  and  surface  micromachining  techniques  are  used  in  the 
fabrication  process.  Burst  pressures  as  high  as  75  psi  have 
been  achieved  for  the  integrated  system  which  are  limited  by 
the  500x500  /Jin2  gas  inlet/outlet  diapliragms. 

Abstract 

Fluid-driven  mechanisms  for  micromechanical  systems  are 
attractive  alternatives  for  microactuators.  Unfortunately,  mi¬ 
cro  flows  with  low  Reynold’s  numbers  are  not  well  under¬ 
stood.  Several  works  have  been  done  to  study  microchan¬ 
nel  flows[l,  2,  3],  but  so  far  there  is  no  experimental  data 
on  pressure  distribution  along  a  microchannel.  This  miss¬ 
ing  information  is  very  important  because  it  can  be  used  to 
develop  flow  models  (since  conventional  fluid  dynamics  the¬ 
ory'  can  not  predict  the  microflow  behavior[3j).  To  provide 
this  information,  an  integrated  microchannel/sensor  system 
for  direct  pressure  measurement,  is  presented  here. 

Fig.  1  shows  a  schematic  of  the  designed  microchan- 
nel/pressure  sensor  system.  The  size  of  the  sensor  diaphragm, 
250x250  pm2  ,  is  selected  to  provide  sufficient  pressure  sensi¬ 
tivity  to  detect  a  pressure  drop  between  any  two  neighboring 
sensors  which  are  500  /im  apart  along  a  microchannel.  A 
small  bridging  channel  between  the  pressure  sensor  and  the 
microchannel  is  also  designed  to  minimize  the  sensor’s  inter¬ 
ference  with  the  flow.  Piezoresistors  on  the  pressure  sensor 
are  connected  to  form  a  Wheatstone  bridge  for  temperature 
compensation  (Fig.  3).  In  addition  to  the  structures  shown 
in  Fig.  1,  venturi  channels  with  pressure  sensors  (not  shown) 
are  also  fabricated  for  velocity-pressure  calibration. 

Fig.  2  shows  the  major  fabrication  steps.  First,  1000 
A  LPCVD  nitride  is  deposited  and  patterned  to  define  di¬ 
aphragms  and  microchannels.  Exposed  silicon  substrate  is 
etched  down  6000  A  and  a  1.0  /mi  thick  oxide  layer  is  then 
grown  by  wet  oxidation.  Next,  a  3000  A  LPCVD  PSG  sacrifi¬ 
cial  layer  is  patterned  (Fig.  2a).  A  1.2  /im  thick  LPCVD  low- 
stress  nitride  (Si,Nj,)  layer  is  deposited  as  the  pressure  sensor 
diapliragm.  A  boron-doped  polysilicon  layer  2500  A  is  then 
deposited  and  patterned  to  form  piezoresisters  for  pressure 
sensors  (Fig.  2b).  After  a  2000  A  nitride  layer  deposition,  gas 
inlet/outlet  holes  are  patterned  on  the  back  of  the  wafer  and 
etched  to  the  front  inlet/outlet  diapliragms  in  EDP  (8  hours). 

Front  holes  for  sacrificial  laver  et.rhine  are  then  onened  by  an 
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SFe  plasma.  High  concentrated  HF  is  then  used  to  etch 
PSG/SiC>2  sacrificial  layer  to  form  micro  channels  and  pres¬ 
sure  sensors  (Fig.  2c).  An  accurate  time-controlled  etch  is 
crucial  in  this  process  [4].  To  seal  the  front  etching  holes,  1.0 
/im  thick  PECDV  nitride  is  deposited.  A  metallization  step 
follows  to  finish  the  process. 

Fig.  5  shows  the  integrated  microchannel/pressure  sensor 
system.  This  system  has  two  channels  (the  top  one  is  5  //m 
wide,  the  bottom  one  is  40  /im  wide)  with  four  sensors  dis¬ 
tributed  along  each  of  them.  The  experimental  setup  of  the 
system  is  shown  in  Fig.  6.  The  mechanical  strength  of  the 
low-stress  nitride  diapliragm  determines  the  maximum  pres¬ 
sure  that  can  be  loaded  to  the  system.  Our  test  shows  that 
the  gas  inlet/outlet  diaphragms,  500x500  pm2  ,  are  the  weak¬ 
est  parts  of  the  system,  but.  they  survive  static  pressures  as 
high  as  75  psi.  We  notice  no  plastic  deformation  of  the  whole, 
system  after  many  pressure  cycles  (65  psi).  A  photograph  of 
the  pressure  sensor  is  shown  in  Fig.  3.  The  resistance  of  each 
wheatstone  bridge  arm  is  18  k fi  .  In  Fig.  4a,  an  SEM  photo¬ 
graph  shows  a  pressure  sensor,  the  bridging  channel,  and  the 
flow  channel.  This  sensor  displays  a  linear  range  from  -15  psi 
to  15  psi.  It  saturates  at  20  psi,  where  the  diaphragm  touches 
the  bottom  of  the  sensor  chamber  with  the  pressure  applied 
from  outside.  Fig.  4b  gives  a  sensitivity  of  120  /iV/V/psi 
which  agrees  well  with  our  design  and  is  enough  to  detect  a 
pressure  drop  between  any  two  neighboring  sensors  500  //m 
apart.  Extensive  measurements  of  the  pressure  drop  along 
microchannels  are  under  way. 
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Fig.  2.  Major  fabrication  steps  of  the  integrated  microchan- 
nel /pressure  sensor  system.  Not  shown  are  the  gas  in¬ 
let  /out  let  holes  etched  from  the  back  of  the  silicon  wafer  for 
fluid  connections  (Fig.  5). 


Fig.  4.  (a)  SEM  photograph  of  a  pressure  sensor,  (b) 

Voltage  output  vs.  pressure  characteristics.  The  bias  voltage 
of  the  Wheatstone  bridge  is  4  V.  The  signal  saturates  when 
the  diapliragm  touches  the  chamber  bottom  at  about  20  psi. 
The  linear  part,  gives  a  sensitivity  of  120  pV/V/psi. 


Fig.  5.  Photograph  of  a  finished  microchannel /pressure  sen¬ 
sor  chip,  ready  for  connections  to  external  gas  sources.  The 
chip  has  two  flow  channels;  the  bottom  is  40  pm  wide  and  the 
top  one  is  5  pm  wide;  Four  pressure  sensors  (250  x  250  pm2  ) 
are  connected  to  each  microchannel  through  a  small  (5x0.3 
pm2  )  bridging  channel  to  monitor  the  pressure  distribution 
along  the  channel. 


Fig.  3.  Photograph  of  a  fabricated  pressure  sensor.  Piezore¬ 
sistors  are  made  of  boron  doped  polysilicon.  Each  arm  resis¬ 
tance  of  the  wheatstone  bridge  is  18  kfl  . 


Fig.  6.  Schematic  experimental  setup  of  the  system.  The 
channel/sensor  chip  is  packaged  in  a  DIP.  The  fluid  connec¬ 
tions  are  through  machined  holes  on  the  bottom  of  the  DIP. 
Electrical  connections  are  through  bonding  wires. 
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Summary  —  A  surface  micromachined  pressure  sensor  using  a  silicon  nitride  diaphragm  and  an  oxide/ PSG  composite 
sacrificial  layer  is  presented(Figs.  1  and  2).  Special  emphases  are  important  micromachining  issues  that  produce  a  high 
diaphragm  yield.  These  issues  include  minimizing  free-etching  time,  avoiding  bubble  formation,  and  designing  proper  etching 
geometry.  Accordingly,  a  diaphragm  yield  of  almost  100%  is  obtained  with  membrane  sizes  smaller  than  200x200  fim7. 

Abstract —  Surface-micromachined  pressure  sensors  are  attractive  for  many  reasons[l].  Guckel  reported  a  polysilicon 
diaphragm  pressure  sensor  using  thermal  oxide  as  a  sacrificial  layer[2].  However,  no  information  on  yield  issues  was  reported. 
Silicon-nitride  pressure  sensors  have  also  been  reported[l,3]  using  silicon  and  polysilicon  as  the  sacrificial  layers,  respectively. 
It  is  pointed  out  that  hydrogen  bubbles  produced  from  silicon  etching  using  KOH  can  cause  membrane  fracture,  hence  reduce 
the  yield.  Here,  we  describe  a  surface-micromachined  silicon-nitride  pressure  sensor  with  a  new  design  that  leads  to  a  high 
diaphragm  yield  of  almost  100%. 


In  our  design,  both  square  and  circular  pressure  sensors  are  available.  Their  sizes  are  150,  200,  250  and  500  fim  on  the 
side  or  in  diameter.  Figure  3  shows  major  steps  of  our  process.  First,  cavity  windows  are  opened  on  silicon  nitride  to  expose 
silicon  substrate.  An  isotropic  silicon  etchant  then  creates  660  nm  step  into  the  substrate.  A  recessed  LOCOS  process  follows 
to  planarize  the  surface  and  form  a  1.3  fim  thick  thermal  oxide  island.  Next,  a  200  nm  thick  PSG  layer  is  deposited  and 
patterned  to  complete  the  composite  oxide/PSG  sacrificial  layer.  This  step  is  followed  by  a  800  nm  LPCVD  nitride  deposition 
to  provide  membrane  material.  Polysilicon  is  then  deposited,  doped,  patterned,  and  passivated  to  form  piezoresistors.  Next, 
an  SF6  plasma  etches  through  the  nitride  layers  to  open  etching  holes.  The  composite  sacrificial  layer  is  subsequently  removed 
using  concentrated  HF  (48%).  This  etching  process  is  specially  video  recorded  for  yield  study.  PECVD  nitride  deposition 
then  seals  the  cavity  at  a  pressure  of  200  mTorr.  Finally,  metalization  follows  and  the  whole  devices  except  the  diaphragm 
pads  is  passivated  using  PECVD  nitride. 


^bonding 


In  our  process,  silicon  nitride  is  chosen  for  the  diaphragm  material  because  of  its  good  mechanical  properties.  However, 
nitride  membranes  may  still  be  damaged  during  fabrication  process.  It  is  found  that  there  are  three  major  causes  of  damage. 
First,  prolonged  HF  etching  attacks  silicon  nitride,  although  at  a  slow  rate.  Secondly,  the  geometry  of  the  etching  holes 
creates  special  free-standing  patterns  of  the  diaphragm  during  etching  process,  which  may  cause  stress  concentration  and 
diaphragm  fracture(Fig.4).  Thirdly,  the  48%  HF  etchant  is  found  to  generate  bubbles  after  the  sacrificial  layer  is  completely 
consumed,  which  may  break  or  lift  off  the  diaphragm.  In  conclusion,  our  studies  show  that  the  solutions  to  these  problems 
are  to  minimize  the  free-etching  time,  to  design  proper  etching  patterns,  and  to  avoid  over  etch.  In  order  to  minimize  the 
etching  time,  an  oxide/PSG  composite  sacrificial  layer  is  used  for  two  major  advantages.  First,  thermal  oxide  allows  a  recessed 
LOCOS  process  to  planarize  the  surface.  Secondly,  the  etching  time  of  the  composite  layer  is  minimized  because  of  the  fast 
etching  rate  of  PSG  in  HF.  For  example,  a  200x200  fim7  diaphragm  cavity(totally  2  fi m  deep)  can  be  freed  in  35  minutes. 
Geometry  of  the  etching  patterns  is  also  an  important  issue.  Sizes,  locations  and  number  of  the  etching  channels  determine 
the  diaphragm  free-etching  time  and  stress  distribution.  Different  etching  channel  patterns  are  included  in  our  design.  It 
is  found  from  in  situ  video  monitoring  that  diaphragm  fracture  happens  at  high  stress  areas  such  as  the  edges  of  etching 
channels  and  the  etching  patterns  during  the  freeing  process  (Figure  4).  Finally,  the  bubble  formation  during  sacrificial  layer 
etching  is  detrimental  but  avoidable  using  a  precisely  timed  etching.  This  suggests  that  uniform  diaphragm  sizes  over  the 
whole  wafer  is  desirable  because  bubble  formation  happens  mostly  after  the  diaphragm  is  completely  freed.  The  diaphragm 
size  effects  on  the  yield  has  also  been  studied.  It  is  found  that  a  yield  of  100%  is  obtainable  for  diaphragm  sizes  smaller  than 
250x250  fim7.  Larger  diaphragms  have  a  samller  yield.  For  example,  the  yield  for  500x500  fim7  diaphragms  is  only  60%. 
However,  the  size  effects  on  the  yield  depend  heavily  on  the  process  and  nitride  quality. 


Finally,  finished  pressure  sensors  are  also  tested.  Typically,  the  output  response  is  linear  within  a  range  from  14  to 
30  psi,  with  a  pressure  sensitivity  of  0.27  mV/V/psi(Fig.  5)  and  a  temperature  sensitivity  of  6/iV/V/°C.  The  polysilicon 
piezoresistors  are  heavily  doped  (>  1020cm-3)  to  have  a  small  temperature  coefficient  of  resistance,  0.02%  /  °C  from  -50  to 
mo  °C  (Fig.  6)  and  a  gauge  factor  of  8.  % 
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Figure  1:  A  square  pressure  sensor.  The  interference 
pattern  is  caused  by  sealing  of  the  diaphragm  at  200 
mTorr. 


Figure  2:  A  circular  pressure  sensor  with  etching 
holes  left  open  on  purpose  for  temperature  sensitivity 
study.  Interference  pattern  therefore  disappears. 
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Figure  4:  A  circular  diaphragm  with  specially  de¬ 
signed  etching  holes  for  etching  pattern  study.  The 
etching  process  was  stopped  half-way  through.  The 
diamond-shaped  area  at  the  center  of  the  diaphragm 
is  the  remaining  axide/PSG  layer.  Sharp  comers  de¬ 
veloped  during  the  etching  process  to  cause  stress 
concentration  and,  sometimes,  diaphragm  fracture. 
The  two  hot-dog  shaped  areas  are  due  to  diaphragm 
warpage. 


1 

\ 

l 


Figure  5:  Pressure  sensor  response.  A  linear  region 
from  14  to  30  psi  is  obtained.  The  pressure  sensitivity 
is  found  to  be  0.127  mV/V/psi.  The  saturation  over 
30  psi  is  because  the  diaphragm  touches  the  bottom 
of  the  cavity. 


Figure  3:  Major  fabrication  steps  for  the  pressure 
sensor,  (a)  a  composite  oxide/PSG  sacrificial  layer 
covered  by  diaphragm  nitride,  (b)  after  diaphragm 
freeing  etching  and  sealing,  (c)  the  complete  struc¬ 
ture  after  metalization. 


Figure  6:  Temperature  sensitivity  characteristic  of 
the  output  voltage  for  an  un-sealed  pressure  sensor. 
The  piezoresistor  Wheatstone  bridge  is  biased  at  4V, 
and  the  reference  voltage  is  0.221  V  at  30°C  due  to 
mismatch  of  the  piezoresistors.  An  averaged  output 
temperature  sensitivity  of  6/iV/V/°C  is  obtained. 
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Silicon  Micromachined  Microstructures  for 
Super-Compact  Magnetic  Recording  Rigid  Disk  Drives 

D.K.  Mill*  P.C.  Chan,1  S.  Tatic-Ludc,1  and  Y.C.  Tais 


Summary 

This  paper  focuses  on  the  design  and  fabrication 
of  millimeter-scale  silicon  microstructures  for  super¬ 
compact  magnetic  recording  rigid  disk  drives  (Fig¬ 
ures  2  to  4).  The  fabrication  processes  consist  of  a 
combination  of  RIB  and  EDP  etching  on  epitaxial 
wafers  with  heavily  boron-doped  buried  layers.  Both 
experimental  data  (Figure  5)  and  finite-element  com¬ 
putations  (Figure  6)  are  presented. 

Abstract 

A  new  generation  of  super-compact  (<  1. 3-inch  form 
factor)  magnetic  recording  rigid  disk  drives  is  being 
introduced  to  accommodate  the  storage  and  process¬ 
ing  requirements  of  real-time,  high-resolution  digi¬ 
tized  graphical  images  in  extremely  challenging  com¬ 
puting  environments  such  as  those  in  multimedia 
portable  computers,  interactive  televisions,  digital 
cameras,  and  handheld  consumer-oriented  informa¬ 
tion  and  communication  systems. 

The  activities  described  here  are  part  of  a  larger 
industry-supported,  jointed  UCLA /CalTech  research 
program  to  design  and  fabricate  a  variety  of  sili¬ 
con  microstructures  and  microactuators  for  super¬ 
compact  data  storage  devices.  Specifically,  our  focus 
here  is  a  micro-gimbal/suspension  using  silicon  as  the 
substrate  material. 

The  gimbal /suspension  assembly  (Figure  1)  is  a 
very  critical  component  for  high-performance  mag¬ 
netic  disk  drives  and  is  designed  to  serve  some  very 
important  functions.  Topically,  one  thinks  of  sili¬ 
con  as  a  brittle  material  and  would  not  normally 
use  it  to  construct  flexible  structures.  However,  it 
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is  well  known  that  the  ultimate  strength  of  silicon 
microstructures  is  orders  of  magnitude  higher  than 
that  in  bulk  form.  Therefore,  when  properly  micro- 
machined,  the  silicon  suspension  can  be  made  very 
flexible  but  with  a  very  high  load-carrying  capacity. 

Fabrications  of  our  micro-gimbals/suspensions  are 
performed  on  4-inch  epi-wafers.  The  composite  epi¬ 
taxial  layers  on  these  wafers  consist  of  a  16  /zm  lightly- 
doped  layer  on  top  of  a  4  /zm  heavily  boron-doped 
buried  layer.  This  buried  boron  layer  is  used  as  an 
etch-stop  for  subsequent  anisotropic  etching. 

A  500  A  chromium  is  first  evaporated  on  the  epi- 
wafers,  which  is  then  patterned  using  photolithogra¬ 
phy,  outlining  the  contours  of  the  microstructures. 
RIE  then  follows,  etching  through  both  epitaxial  lay¬ 
ers.  After  removal  of  any  residual  chromium,  wafers 
are  either  wet-oxidized  to  a  1  /zm  thermal  oxide  or 
deposited  with  a  4  /zm  PECVD  oxynitride  to  form 
protective  layer  on  all  surfaces.  Second  photolithog¬ 
raphy  step  is  performed  on  the  backside  of  the  epi- 
wafers  enabling  the  opening  of  windows  in  the  oxide 
(or  oxynitride)  layer.  EDP  etching  is  then  performed 
to  form  cavities  under  the  microstructures .  Finally, 
oxide/nitride  is  removed  to  release  the  desired  free¬ 
standing  microstructures. 

Figures  2  to  4  present  the  SEM  pictures  of  a  num¬ 
ber  of  prototype  micro-gimbals,  measuring  typically 
2.5  by  2.5  mm.  Specifically,  Figure  2  shows  a  scaled- 
down  gimbal  design  resembling  that  of  the  industry- 
standard  type- 16  suspension.  Figure  3  is  a  simpler  de¬ 
sign  consisting  of  only  straight  cantilever  beams  and 
has  a  form  factor  consistent  with  the  so-called  35% 
sliders.  Finally,  Figure  4  shows  some  test  structures ' 
for  the  purpose  of  determining  the  micromechanic&l 
properties  of  silicon. 

Component-level  performance  characteristics  of 
these  microstmctures  are  investigated,  using  both  a 
ultra-sensitive  load-cell  for  static  deflection  and  frac¬ 
ture  measurements  as  well  as  a  nonintrusive  laser 
Doppler  vibrometer  for  dynamic  analysis  (Figure  5). 
In  addition,  finite-element  calculations  are  performed 
(Figure  6),  allowing  comparison  and  design. 
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Figure  1:  Detailed  drawings  of  the  conventional  3370- 
type  read/write  suspensions  with  a  screw  and  a  swage 
mount. 


Figure  2:  SEM  micrograph  of  a  type- 16-like  silicon 
micromachined  gimbal  (top  view). 


Figure  3:  SEM  micrograph  of  a  simple-frame  silicon 
micromachined  gimbal  (bottom  view). 


Figure  4:  SEM  micrograph  of  some  simple  test  mi- 
crostructures  (top  view). 


Figure  5:  Frequency  response  of  the  simple  test  struc¬ 
ture  using  a  nonintrusive  laser  Doppler  Vibrometer. 


Figure  6:  Finite  element  analysis  of  the  simple- frame 
silicon  gimbals. 
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(3.  Fabrication  Technologies) 

Mass-Producible  All- Silicon  Atomic  Force  Microscope  Probes 
Chang  Liu  and  Yu-Chong  Tai 

Electrical  Engineering  116-81,  Caltech,  Pasadena,  CA911S5,  USA 
Ronald  Gamble 

Topometries  Corp.,  Santa  Clara,  CA95054-1524,  USA 

Summary  —  Mass-producible  all-silicon  Atomic  Force  Microscope  (AFM)  probes  (Figures  1  and  2)  are 
presented.  A  new  approach  of  using  heavily  boron-doped  buried  layer  as  the  beam  material  to  provide  robust 
fabrication  process  is  adopted.  Mechanical  characteristics  of  the  probes  have  been  measured  by  a  Doppler  laser 
vibrometer.  Good  AFM  images  (Fig.  4)  are  obtained  using  these  probes. 

Abstract  —  Since  its  invention,  the  Atomic  Force  Microscope  (AFM)  has  become  a  powerful  tool  for  surface 
study  with  atomic  resolution1.  However,  in  order  to  maintain  its  high  performance,  it  is  mandatory  to  have 
high  quality  probes  with  reliable  mechanical  characteristics  and  good  tip  sharpness.  Here,  we  then  present  an 
all-silicon  probe  which  meets  these  requirements.  This  probe  has  three  parts:  a  tip,  a  beam  and  a  handle.  The 
probe  is  completely  made  of  single  crystalline  silicon  in  order  to  have  highly  reliable  mechanical  properties3.  To 
make  the  probe,  we  design  a  process  using  (100)  epitaxial  wafers  that  have  a  10  fim  lightly  doped  layer  on  top  of 
a  4  n m  heavily  boron-doped  buried  layer.  The  lightly  doped  layer  is  used  to  form  the  sharp  tips,  and  the  heavily 
doped  buried  layer  is  to  form  the  beam.  This  approach  of  using  a  buried  layer  to  make  the  beam  allows  good 
control  of  the  beam  thickness  since  EDP  etching  automatically  stops  at  this  layer.  To  our  knowledge,  this  is  the 
first  time  this  approach  has  been  reported. 

Major  steps  of  our  process  is  described  here.  The  fabrication  starts  with  forming  silicon  pyramids  in  the 
lightly  doped  epitaxial  layer  by  EDP  etching  (Figure  6a).  A  l-/im-thick  low-temperature-oxide  (LTO)  masking 
layer  is  then  deposited  on  the  wafer  and  patterned  to  define  the  cantilever  beams.  The  beams  are  etched  by 
SFg  plasma.  An  LPCVD  low-stress  nitride  layer  (200  nm  thick)  is  deposited,  followed  by  beam-freeing  EDP 
etch  from  the  back  of  the  wafer.  Next,  the  LPCVD  nitride  is  stripped  in  phosphoric  acid  to  free  the  beam.  Tip 
sharpening  then  finishes  the  process  by  growing  wet  oxide  on  the  tip  and  removing  it(Figure  6b).  The  finished 
tips  typically  have  a  radius  of  less  than  100  nm  and  an  aspect  ratio  of  1.  Figures  1  and  2  show  a  probe  with  a 
sharp  tip,  a  V-shaped  beam  and  a  handle.  The  beam  is  200  /im  long,  20  pm  wide,  and  2  pm  thick. 

Mechanical  characteristics  of  the  probes  have  been  measured  using  a  vibrometer,  which  consists  of  a  Doppler 
laser  interferometer.  For  example,  the  resonant  frequency  of  the  probe  shown  in  Figure  1  is  42  KHz ,  which  is  in 
good  agreement  with  finite  element  simulations.  Quality  factor,  Q,  of  the  probe  is  also  found  to  be  80  in  the  air 
using  both  band-width  and  time-domain  step  response  methods  (Figure  3).  Aging  of  the  boron-doped  beams 
has  also  been  studied  by  exciting  the  beam  at  its  resonant  frequency  for  as  long  as  10  hours  (about  l.SxlO9 
cycles),  and  no  measurable  change  of  the  resonant  frequency  and  Q  was  found. 

The  probes  have  also  been  tested  in  AFM  machines  (Topometries  TMX-2000)  under  both  contact  and  non- 
contact  modes.  Figure  4  shows  an  AFM  image  of  a  compact  disc  surface  using  our  probes.  As  shown,  the 
pmrsized  pits  are  clearly  seen.  Finally,  a  M Superprobe”  process  can  be  added  to  our  probes  so  that  even  better 
tip  sharpness  and  aspect  ratio  can  be  obtained,  if  necessary.  Figure  5  demonstrates  such  a  superprobe.  The 
radius  of  the  superprobe  is  about  10  nm  and  the  aspect  ratio  is  more  than  10. 
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Figure  1:  SEM  picture  of  a  V-shaped  beam.  At  the 
end  of  the  beam  is  an  oxidation-sharpened  tip. 


Figure  2:  SEM  picture  of  a  complete  probe  which 
consists  of  a  tip,  a  beam  and  a  handle.  The  handle  is 
designed  for  easy  mounting  of  the  probe  to  AFMs. 


Figure  4:  An  AFM  image  of  a  compact  disc  surface 
obtained  using  the  probe  shown  in  Fig. 2.  The  CD 
track  pitch  is  about  2.5  /xm  and  the  pits  are  ljxm  in 
diameter. 
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Figure  5:  SEM  picture  of  a  superprobe  grown  on  a 
silicon  tip. 
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A  variety  of  flow  structures  have  been  suggested  as  important  in  the  development 
and  regeneration  of  turbulent  boundary  layers.  The  majority  of  these  suggested  flow 
structures  are  associated  with  some  form  of  vortex  structure,  from  streamwise  and 
transverse  vortices  to  multiple  hairpin  vortices.  However,  most  of  the  flow  structures 
and  their  hypothesized  behavior  are  based  on  inferences  from  measurements  or  simula¬ 
tions  of  fully-turbulent  flows,  often  revolving  around  a  discussion  of  kinematics,  but 
with  little  dynamic  support.  A  particularly  important  process  involved  in  the  dynamics 
of  turbulence  is  that  of  vortex  interaction  with  surfaces,  a  process  which  provides  a 
sustaining  mechanism  for  turbulence  through  the  regeneration  of  new  vortices.  The 
present  talk  will  summarize  the  findings  of  a  series  of  ''kernel"  experiments  which 
have  examined  the  basic  kinematic  and  dynamic  behavior  of  various  vortex  phenomena 
believed  to  be  relevant  to  wall  turbulence,  focussing  on  how  vortices  interact  with  sur¬ 
faces,  and  how  these  interactions  lead  to  generation  of  "new"  vorticity  and  vortices. 
Using  both  flow  visualization  and  PIV  results,  several  different  vortex  flow  structures 
will  be  discussed,  including  ring  vortices,  streamwise  vortices,  hairpin  vortices,  and 
necklace/horseshoe  vortices  in  junction  regions.  A  conceptual  model  of  riear-wall  tur¬ 
bulence,  based  upon  the  results  of  these  vortex  interaction  studies  and  other  relevant 
work  will  be  discussed. 
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Intoduction 


Turbulence  —  what  is  it? 

•  Measurements  reflect  property  fluctuations 

•  Originally:  Statistical  interpretation 

•  Recently:  Coherent  structure  interpretation 

Difficulty  =>  Establishing  and  Modeling  of 
relevant  flow  structure 

Example:  Flow  visualization  of  TBL  on  flat  plate 

•  H2  bubbles  in  water  flow 


•  Side-view  &  Plan-view 


Suggestion  of  organization  =>  Coherent  Structure 

•  Complicated,  unclear  interactions 

Kucheman  =>  ’Vortices  are  sinews  and  muscles  of  turbulence" 

•  What  part  do  vortices  play? 

•  Where  do  they  orginate? 

•  How  do  they  interact: 

-  with  each  other? 

-  with  the  surface? 
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Vortex  Kernel  Studies 


Kernel  study: 

•  Examination  of  a  simpler  flow  structure  in  order  to  understand 
the  dynamics  and  processes  relevant  to  turbulence 

Examples: 

•  Deformation  of  aline  vortex  in  a  shear  flow 

•  Interaction  of  a  controlled  vortex  with  a  surface 

Objective: 

•  To  understand  the  processes  giving  rise  to  observed/detected 
turbulent  behavior 

To  be  considered  in  this  seminar 

•  Review  of  selected  kernel  studies 

V  Inviscid  3-D  vortex  interactions  in  shear  flow 

V  Vortex-surface  interactions 

V  Hairpin  vortex  evolution/interactions 

•  Video  of  selected  behavior 

•  Wall  turbulence  model  suggested  by  kernel  studies 
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Inviscid  Vortex  Deformation/Interaction 


(Courtesy  of  JJ3A.  Walker  et  aL) 
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Fig.  3.  Primary  structure  of  wallbound 
turbulence  (Presented  at  the  Secound 
Midwestern  Conference  on  Fluid  Me* 
chanics,  March  17,  1952) 


Fig.  4.  Horseshoe  of  large  Rtyno'.di  number 
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Initial  vortex  position 


Evolution  of  a  small  step”  in  a  two-dimensional  vortex  convected  in  a 
uniform  shear  flow:  (a)  top  view,  (b)  side  view,  (c)  end  view;  £  =  0.00267, 
440  time  steps  with  At  =  0.034.  The  vortex  position  is  plotted  every  40 
time  steps.  pc 

TY^ HA  4 


Initial  vortex  position 


Initial  vortex 


Initial 


vortex 


Figure  7.  Evolution  of  a  “step”  in  a  two-dimensional  vortex  convected  in  a  uniform 
shear  flow;  the  vortex  position  is  shown  every  40  time  steps  in  (a)  top  view, 
(b)  side  view,  and  (c)  end  view.  , 


^  <1  4 


Figure  12.  Evolution  of  multiple  hairpin  vortices  in  a  uniform  shear:  (a)  top  view,  (b) 
side  view,  (c)  end  view;  all  three  coordinates  are  plotted  on  the  same  scale. 


II- 


Head  Subsidiary 


Figure  6.  Evolution  of  a  symmetric  hairpin  vortex  in  a  shear  flow. 

(a)  initial  distortion,  (b)  development  of  vortex  legs  and  head, 
(c)  evolution  of  subsidiary  vortices  and  penetration  toward  the 
surface. 
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(b)  Asymmetric 


;ively,  expect  irregularity 


s  of  interaction  similar 


Vortex-Surface  Interaction 
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(a)  The  evolution  of  recirculation  (b)  Eddy  growth 

near  the  surface 


(c)  Focussing  and  compression;  solid  lines  denote  the 
equi-vorticrty  contours  in  the  surface  layer;  arrows 
indicate  compression  of  the  erupting  spike 


) 


Figure  3.  Schematic  diagram  of  the  stages  in  the  separation  of  a 
two-dimensional  surface  layer  (not  to  scale). 
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counter-rotating 


wall 


focusing  of 


wall 


1)  An  isolated  vortex  is  in  the 
vicinity  of  a  solid  wall. 


2)  A  counter-rotating  eddy 
develops,  and  surface  fluid 
lifts  from  the  wall. 


3)  The  combined  influences  of  the 
parent  vortex  and  the  strengthening 
counter-rotating  vortex  have  a 
focusing  effect  on  the  surface  fluid 
lifting  off  the  wall. 


4)  The  narrow  band  of  surface  fluid 
violently  erupts  away  from  the  wall, 
interacting  with  the  outer  flow.  The 
counter-rotating  eddy  is  ejected 
into  the  outer  flow  during  the 
process. 


Schematic  Illustration  of  a  Vortex-Surface  (Inviscid-Viscous)  Interaction. 
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Necklace 

Vortices 


Bluff 

Body 


Processes  in  Unsteady  End-Wall  Flow: 
Turbulent  ADDroach  Flow 


End  Wall 


Impinging  Turbulent 
Boundary  Layer 
Structure 


Vortex  Generation 
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sv  I  (C)  (d) 


(eruption) 


Figure  14.  Four  sequential  particle  image  visualizations  of  surface  interaction  on  symmetry 
plane  of  turbulent  horseshoe  vortex  injunction  region  of  rectangular  body. 
Re0=7OO;  each  picture  taken  0.25  seconds  apart. 

HV=*horseshoe  vortex;  SV=*secondary  vortex. 
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PIV  REALIZATIONS  OF  A 
TURBULENT  HORSESHOE  VORTEX 


STREAMLINES  VORTICITY 


O. 
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Light  Sheet 
Illumination 


(c)  Induced  surface  eruptions 


(d)  Rollover  into  asymmetric  hairpin 


(e)  Lagrangian  computation  of  vortex-induced, 
temporal  development  of  a  materia!  line 
illustrating  development  of  eruptive  spike 


Figure  1.  End-view  of  surface  eruptions  and  hairpin  generation  induced  by 
necklace  vortex  legs.  Visualization  using  hydrogen  bubbles  and 
light-sheet  illumination. 
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Plan-view  photograph  showing  a  streamwlse  necklace  vortex 
leg  (visualized  by  upstream  hydrogen  bubble  wire)  interacting 
with  suilace  fluid  (visualized  by  bubble  wire  near  surface)  to 
create  low-speed  streak  at  surface. 


Figure  1 8  Plan-view  of  necklace  vortex  leg  interaction  to  create  low-sp^ed  streak 
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Low-Speed  Streaks: 
Turbulent  Boundary 

Layer,  Ke6  =1500 


Surface  Eruption  Caused 
by  Streamwise  Laminar, 
Junction  Vortex 


Comparison  of  Plan- View  Hydrogen  Bubble 
Visualization  of  Flow  Near  the  Surface  Beneath: 

A)  a  Turbulent  Boundary  Layer,  and 
B)  the  Streamwise  Portion  of  a  Necklace  Vortex 
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hairpin  streamwise 

vortex  vortex 


a)  End-view,  light-sheet  photograph  showing  the  breakdown  of  an  eruptive 
spire  (low-speed  ridge)  induced  by  a  streamwise  necklace  vortex  leg. 


b)  Schematic  illustrating  the  general  breakdown 
process  for  the  low-speed,  eruptive  ridge. 


Figure  20.  Characteristics  of  the  breakdown  of  the  eruptive  spires  induced  by 
streamwise  necklace  vortex  legs. 
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Hydrogen 


Figure  25.  Schematic  diagram  of  laser  sheet  visualization  technique  employed 
to  visualize  eruptive  events  in  a  turbulent  boundary  layer. 


(a) 


(b) 


Figure  19.  End-view  visualizations  of  surface  eruptions  in  both  a  transitional  and  turbulent 
boundary  layer,  (a)  Thin,  eruptive  spire  adjacent  to  surface  in  a  tripped, 
transitional  boundary  layer  (hydrogen  bubble  visualization),  (b)  Series  of  eruptive 
surface  spires  in  a  turbulent  boundary  layer  (surface  smoke  visualization)  from 
Wallace  et  al.  (1990). 


11-34 


I)  It 


Hairpin  Vortex  Evolution/Interaction 
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Early  Stages  Of  Observation 
Perry,  Lim  and  Teh  (1981) 


LEGS 


PLAN-VIEW 


FLOW  - 


SIDE-VIEW 


(a) 


* 


Head 


Jigure  15.  Illustration  of  experimental  generation  of  a  single  hairpin  vortex  by 

surface  injection  using  system  shown  in  figure  14.  (a)  Dual-view  picture 
of  dye-marked  single  hairpin  vortex,  (b)  Isometric  schematic  of  single 
hairpin  after  generation. 
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Viewing  Direction: 
Figures  17  and  23 


Hydrogen 
Bubble  Probe 


Hydrogen 

Bubble 

Sheet 


Hydrogen  Bubble  Pattern 
Created  by  Single,  Advecting 
in  Vortex 


Hairpin  Vortex 
(dark  portion 
not  visualized 
in  bubble  sheet) 


Illumination  from  Below 


Figure  1 6.  Schematic  diagram  of  technique  employed  to  visualize  flow 
development  associated  with  a  single  hairpin  vortex 
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Dual-Level,  Plan-View  Visualization  of  Hydrogen  Bubble  Patterns 
Generated  by  Developing  Hairpin  Vortex 
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FLOW  ►► 


JHBW  jHBW 


regions. 


x/5=1.0 


Figure  4.16 


x/8=5.0 

(Continued)  for  caption  see  figure  4.16. 


x/8=10 


x/5=30 


FLOW  ►► 


Primary  Hairpin  Vortex 


►Increasing  Time 


e)  t=  1.0  0  t=  1.25  g)  t=  1.50  h)  t=  1.75 

Figure  5.1  Side-view  dye  visualization  of  a  turbulent  spot  generated  from  growth  of  a  single  hairpin  vortex.  Photographs  are 
obtained  with  still  photograplty  at  4  frames  per  second.  Reg  =440,  ReVw=10.5,  x=  60cm. 
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b’)  x=7.0  cm 

y=0.25  cm 

Sp*h  CT.6.) 


b)  x=1.0  cm 

y=0.39  cm 

Vo*-*** 


"  \ 
Figure  5.6 


Effect  of  Reynolds  number  on  the  growth  to  a  turbulent  spot,  as  visualized  by  hydrogen  bubbles,  a)  through  d) 
Re5i*=s440,  a’)  through  d’)  Reg2*=530. 
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FLOW  ►► 


Figure  5.2 


d>  t=0.25  e)  1=0.333  f)* 

Plan-view  visualization  of  a  turbulent  spot  generated  by  a  haiipin  vortex  The 
vortices  within  the  turbulent  spot.  Re8*=530,  ywirc=  0.57cm,  xwire=7cm.  ’ 


t=  0.416 

sequence,  illustrates  fonnation  of  new 


FLOW  ►► 


) 


a’)  x=7.0  cm 

y=0.57  cm 


a)  x=1.0  cm 

y=0.65  cm 


b’)  x=7.0  cm 

y=0.25  cm 


b)  x=1.0  cm 

y=0.39  cm 


Figure  5.3:  Comparison  between  flow  patterns  for  a  single  vortex  and  local 
flow  patterns  observed  in  a  turbulent  spot,  a),  and  b)  single 
HPV,  Re5*=440,  KeVw=10.5,  xwire=  1  cm.  a')  and  b')  trailing 
vortices  in  a  turbulent  spot  Re6*=530,  ReVw=  7.8  ,  2^,=  7  cm. 
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Hairpin  Vortex 
Formation 
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Interactions 
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Dissipation 


Generation 


Model  of  Generalized  Turbulent  Boundary  Layer  Processes 


What  Kernel  Studies  Indicate 


Low-Speed  Streaks 

•  Caused  by  vortex-surface  interaction 

•  Inviscid-viscous  interaction  =>  other  3-D  vortices 

Propagation  of  Turbulence 

•  Due  to  3-D  vortex  deformation  in  velocity  gradient,  and 

•  Surface  eruptions  resulting  in  secondary  vortex  generation 

Growth 

•  Amalgamation  and  coalescence  of  vortex  scales 

•  Energy  input  =>  vortex  stretching  by  velocity  gradient 
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II-2.  The  Interaction  of  Coherent  Vortices 
with  Solid  Surfaces 
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NUWC  Division  Newport,  R.  I. 

SEMINAR  NOTICE 


THE  INTERACTION  OF  COHERENT  VORTICES 
WITH  SOLID  SURFACES 

Professor  Demetri  P.  Telionis 
Virginia  Polytechnic  Institute  and  State  University 


Experimental  and  numerical  methods  have  been  employed  to  study  the  for¬ 
mation  of  coherent  vortices  generated  over  sharp  edges  and  their  subsequent 
interaction  with  blades  further  downstream.  Laser-Doppler  velocimetry  and  flow 
visualization  were  employed  in  a  water  tunnel.  Periodic  flows  were  generated  and 
two-dimensional  fields  were  obtained  by  conditional  averaging.  The  numerical 
methods  were  based  on  discrete  vortex  dynamics.  The  results  indicate  that  the 
vortex  cores  are  created  by  purely  inviscid  mechanisms  and  should  therefore  be 
termed  ’inviscid  cores’  rather  than  ’viscous  cores’.  When  coherent  structures 
encounter  solid  surfaces,  their  distortion  and  further  development  conform  with 
Kelvin-Helmholtz  theorem.  However,  the  large  pressure  gradients  induce  separa¬ 
tion  and  the  creation  of  separated  bubbles  which  drift  together  with  the  disturb¬ 
ing  vortices. 


Thursday ,  27th  May  1993 
Conference  Room ,  Bldg.  679 
Time:  10:30  AM 

POC:  Dr.  Promode  R.  Bandy opadhyay  (Code  8234;  x2588) 
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Figure  2.12  Traversing  mechanism 
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Figure  2.3.6  The  mirror  tower  assembly 
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Figure  3.4.  Velocity  waveforms  at  four  downstream  locations  in  the  wake:  xjct  =  1.85,  1.925,  2.0, 
2.075;  ylca~  0.1;  k  =  2.0 


Chapter  3 


Figure  3.5.  Velocity  vector  Held  constructed  from  a  single  column  of  measurements:  xjcg  *=  2.0, 
k  *=  2.05,  Re  =  14000 
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Figure  3.6.  Velocity  vector  field  in  reference  frame  of  vortex  (from  Fig.  3.6):  xfct  m  2.0,  k  =  2.05, 
Re  =14000 
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Figure  3.8. 


A  comparison  or  grid  and  column  data:  Total  circulation 
:  xfct  -  2.0,  k-  2.05,  Re  =  14000 


Chapter  3 


•  Rolling  up  of  vortex  sheets  is  independent  of  turbulent 
character. 


•  There  are  no  “viscous  cores.”  Cores  are  created  by  convection, 
hence:  “inviscid  cores.” 


•  Many  problems  involving  coherent  vortices  -  laminar  or 
turbulent  -  could  be  solved  with  inviscid  methods. 


•  New  vorticity  transport  process:  vortex  stretching  normal  to 
vortex  lines. 
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Figure  2.2.1  The  pitching  airfoil  and  target  airfoil 
arrangement  in  the  test  section 
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Figure  4.11.  Vorticitv  contour  maps:  (a)  a  =  0°,  Ali  =  0.09c  (b)  a-10°,  Ali-O.lc  (c)  a"*10, 

Ah  =  0.125c 
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Figure  4.11.  continued:  (a)  a  =  0°,  Ah  =  0.09c  (b)  a 
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Figure  4.12.  continued:  x/c  =  0.15:  (a)  a  =  0°,  Ah  =  0.09c  (b)  «=10°,  Ah  =  0.1c  (c)  a*=10“, 

zMi®  0.125c 
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Figure  4.13.  continued:  x/c  =  0.15 
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•  Greatest  effects  observed  for  skimming  vortex  and  loaded 
blade. 


•  Direct  impact  results  in  “disintegration”  of  the  vortex. 


•  Coherence  of  vortex  is  preserved  if  secondary  vortex  is  created. 


•  Pressure  spikes  weaken  as  vortices  drift  downstream  over  the 
blade. 
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Fig.  1.  Profile-plate  arrangement. 
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Vorticity  conlourB  for  the  25.4  min  plate,  Ah  =  10  mm:  Phase,  (j>  -  ip  • 

3G0° .  Thick  and  thin  lines  represent  positive  (counterclockwise)  and 
Fig-  5.  negative  vorticity  respectively.  The  levels  in  positive  and  negative  contour 
lines  correspond  to  dimensionless  values  of  2,0,10,15,20,30  and  50  starting 
from  the  open  domain.  Spatial  uncertainty  estimate  of  contours:  ±  0.01  C  in 
both  directions.  tt  mo 


) 


•  Smaller  plates  are  more  effective  vortex  “choppers.” 


•  Secondary  vortex  could  be  more  powerful  than  primary 


•  Plates  with  rounded  rather  than  sharp  edges  are  more  effective. 
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NUWC  Division  Newport,  R.  I. 

SEMINAR  NOTICE 


ilIGH  RESOLUTION  VORTEX  SIMULATIONS  OF 

BLUFF  BODIES 

Professor  A.  Leonard 
California  Institute  of  Technology 


The  main  theme  of  this  ongoing  research  is  to  study  the  vortex  dynamics  of  bluff  body  flows  and 
the  relation  of  these  dynamics  to  the  unsteady  forces  experienced  by  the  body.  In  addition, 
means  of  controlling  these  forces  and  dynamics  are  being  explored.  Our  computational  approach 
is  to  use  discrete  representations  of  the  vorticity  transport  equation,  a  form  of  the  Navier-Stokes 
equations  for  incompressible  flows.  In  particular,  computational  elements  of  vector-valued 
Vorticity  are  convected  with  the  local  fluid  velocity  and  the  vorticity  vector  associated  with  an 
element  is  continually  strained  by  the  local  velocity  gradient.  A  mix  of  improvements  to  the  usual 
vortex  method  were  implemented  for  these  studies,  including  a  fast  vortex  algorithm  and 
improved  treatment  of  viscous  effects.  In  two-dimensions,  unsteady,  separated  flows  past  a  flat 
plate,  nominally  at  90  degrees  angle  of  attack,  and  past  a  cylinder  are  investigated,  using  a 
Cray-YMP.  Comparison  with  tow-tank  experiments  is  made.  For  the  three-dimensional  work  we 
use  the  Intel  Delta  parallel  computer.  Extension  of  the  the  viscous  technique  to  three-dimensional 
flows  in  the  presence  of  arbitrarily-shaped  solid  bodies  is  complete  and  implementation  of  this 
extension  with  a  fast  three-dimensional  vortex  algorithm  on  the  Delta  to  study,  for  example,  the 
unsteady  pitching  and  yawing  of  a  three-dimensional  bluff  body  is  nearly  complete. 
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High  Resolution  Vortex  Simulations 
of  Bluff  Body  Flows 

A.  Leonard 

Graduate  Aeronautical  Laboratories 
California  Institute  of  Technology,  Pasadena,  CA  91125 


Physical  Theory 

Navier-Stokes  Equations,  Vorticity  Transport  Equation,  Vortex  Dynamics,  Surface  Generation  of 
Vorticity,  Viscous  Transfer,  Body  Forces  and  Moments 

Numerical  Techniques 

Vortex  Particles,  Remeshing,  Panel  Method,  Particle-Particle  Circulation  Exchange, 
Panel-Particle  Tranfer 

'Parallel  Implementation 

Fast  Algorithm,  Error  Control,  Data  Structure,  Load  Balancing 
Turbulence 

Subgrid-scale  Modeling,  Transition  to  Turbulence 
Applications 

Cylinder,  Sphere,  Prolate  Spheriod, . 
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Vorticity  shedding  model 


i 

i 

i 


11-122 


Tin*  :  0.3 
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NON-OSCILLATING  PLATE 
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Comparison  of  Cp  to  Experiments 


Verification 


) 

•  Starting  vortices  -  Analytical  solution 

(Rott  1956,  Puffin  1978) 

•  Recirculating  bubble  -  Experiments 

(Taneda  and  Honji  1971) 

'V 

•  Ramp-start  flow  -  Current  experiments 
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Comparison  of  circulation  flux  to  analytical  solution 
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(a)  Yorticsty 

(CasQ^utations) 


(b)  Streamlines 
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(c)  Streaklines 
(Experiments) 
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(a)  Vorticity 
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Figure  8.  Shadow-photographs  of  the  interaction  of  an  =  1.25  shock  wave  moving  from  right 
to  left  over  a  spherical  helium  volume  (4.5  cm  diameter),  (a)  t  —  20  ps,  (6)  82  ps,  (c)  145  ps,  (d) 

223  ps,  (e)  350  ps,  (/)  600  ps,  (g)  1504  ps. 

a  shock-on-shock  interaction,  which  forms  a  quadruple  shock  intersection.  At  72  ps 
(figure  7 d),  the  secondary  transmitted  waves  are  seen  outside  the  volume  tangentially 
connected  to  the  (primary)  transmitted  wave  and  connected  at  the  interface  to  the 
internal  reflected  wave,  now  approaching  its  focus.  A  wave  reflected  from  the 
transversal  beam  is  seen  at  bottom  right.  At  82  ps  (figure  7c),  the  internal  reflected 
wave  has  become  a  diverging  wave,  having  passed  through  its  caustic.  Outside,  the 
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Kenneth  S.  Breuer 

Department  of  Aeronautics  and  Astronautics,  Room  33-214 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
(617)  253-1980 
breuer@helmholtz.mit.edu 


Employment 

•  6/91  -  Present.  Charles  Stark  Draper  Assistant  Professor  of  Aeronautics  and  Astronautics  Mas¬ 
sachusetts  Institute  of  Technology,  Cambridge  MA.  Responsible  for  undergraduate  and  graduate  level 
teaching,  advising  and  research.  Current  research  projects  include:  Active  control  of  laminar-turbulent 
transition  using  flexible  structures,  active  control  of  turbulent  boundary  layer  structure.  Design  and 
testing  of  micro-machined  probes  for  fluid  mechanics,  investigation  of  physics  of  ice-accretion  on  air¬ 
foils,  effects  of  roughness  on  laminar-turbulent  transition,  Analysis  and  numerical  simulation  of  bypass 
mechanisms  for  transition  in  boundary  layer  and  pipe  flows. 

•  9/90  -  6/91.  Assistant  Professor  of  Aeronautics  and  Astronautics  Massachusetts  Institute  of  Technol¬ 
ogy,  Cambridge  MA. 

•  10/88  -  8/90.  Post  doctoral  research  associate  Center  for  Fluid  Mechanics,  Turbulence  and  Compu¬ 
tation,  Brown  University.  Research  projects  included:  Error  analysis  in  Chebyshev  computations. 
Computation  of  generalized  eigenfunctions  for  linear  and  nonlinear  problems.  Analysis  of  probability 
distributions  in  turbulent  shear  flows. 

Education 

•  6/88  PhD,  Massachusetts  Institute  of  Technology,  Department  of  Aeronautics  and  Astronautics.  Thesis 
title:  “The  development  of  a  localized  disturbance  in  a  boundary  layer.”  Thesis  committee:  M.T. 
Landahl,  J.H.  Haritonidis  &  S.E.  Widnall. 

•  9/84  MS,  Massachusetts  Institute  of  Technology,  Department  of  Aeronautics  and  Astronautics.  Thesis 
title:  “An  Experimental  Investigation  of  Transitional' Pipe  Flow.” 

•  6/82  ScB,  Brown  University,  Division  of  Engineering.  Graduated  Magna  Cum  Laude  with  Honors. 
Concentration  in  Fluids  and  Thermal  Sciences. 

Current  research  interests 

Experimental,  theoretical  and  numerical  fluid  mechanics,  with  a  special  interest  in  stability,  transitional 
flows  and  coherent  structures  in  fully  turbulent  flows.  Active  control  of  transitional  and  turbulent  boundary 
layers;  turbulent  noise  and  drag  reduction.  Computational  Fluid  Dynamics.  Experimental  techniques  and 
micro-sensor  development. 


Miscellaneous 

•  Invited  researcher  at  NASA/Stanford  Center  for  Turbulence  Research,  Royal  Institute  of  Technology, 
Stockholm,  NASA  Langley  Research  Center. 

•  Reviewer  for  Journal  of  Fluid  Mechanics,  Physics  of  Fluids,  Sensors  and  Actuators,  AJAA  Journal , 
Science 

•  Member  of  APS,  AIAA. 

•  Citizenship:  USA. 
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DR.  MOHAMED  GAD-EL-HAK 
SHORT  BIOGRAPHICAL  SKETCH 


Professor  Gad-el-Hak  received  his  B.Sc.  (summa  cum  laude)  in  mechanical 
engineering  from  Ain  Shams  University  in  1966  and  his  Ph.D.  in  fluid  mechanics  from  the 
Johns  Hopkins  University  in  1973,  where  he  worked  with  Professor  Stanley  Corrsin.  He 
has' since  taught  and  conducted  research  at  the  University  of  Virginia,  University  of 
Southern  California  and  Institut  National  Polytechnique  de  Grenoble,  and  has  lectured 
extensively  at  seminars  in  the  United  States  and  overseas.  Dr.  Gad-el-Hak  is  currently 
Professor  of  Aerospace  and  Mechanical  Engineering  at  the  University  of  Notre  Dame. 
Prior  to  that,  he  was  a  Senior  Research  Scientist  and  Program  Manager  at  Flow  Research 
Company  in  Seattle,  Washington,  where  he  has  managed  a  variety  of  aerodynamic  and 
hydrodynamic  research  projects.  Professor  Gad-el-Hak  has  developed  several  novel 
measurement  techniques  for  turbulent  flows,  and  he  holds  two  patents:  one  for  a  drag- 
reducing  method  for  airplanes  and  underwater  vehicles  and  the  other  for  a  lift-control 
device  for  delta  wings.  Dr.  Gad-el-Hak  has  published  over  170  articles  in  the  basic  and 
applied  research  areas  of  isotropic  turbulence,  boundary  layer  flows,  stratified  flows, 
compliant  coatings,  unsteady  aerodynamics;  biological  flows,  non-Newtonian  fluids,  and 
flow  control.  He  is  the  editor  of  two  volumes  in  Springer- Verlag's  Lecture  Notes  in 
Engineering  entitled  " Frontiers  in  Experimental  Fluid  Mechanics and  "Advances  in  Fluid 
Mechanics  Measurements  "  Professor  Gad-el-Hak  is  a  member  of  the  American  Chemical 
Society,  an  associate  fellow  of  the  American  Institute  of  Aeronautics  and  Astronautics,  a 
life  member  of  the  American  Physical  Society,  and  a  member  of  the  American  Society  of 
Mechanical  Engineers.  From  1988  to  1991,  Dr.  Gad-el-Hak  served  as  Associate  Editor  for 
the  AIAA  Journal .  He  is  currently  serving  as  Associate  Editor  for  the  Applied  Mechanics 
Reviews  as  well  as  a  Contributing  Editor  for  Springer-Verlag's  Lecture  Notes  in 
Engineering.  From  1988  to  1990,  Dr.  Gad-el-Hak  served  as  a  member  of  the  Executive 
Committee  of  the  Faculty  Senate  at  the  University  of  Notre  Dame,  Chairman  of  the  Benefits 
Committee,  as  well  as  Chairman  of  the  University  Governance  Committee.  He  served  as  a 
consultant  to  the  U.S.  Government,  the  Government  of  Egypt,  the  United  Nations,  and 
several  industrial  organizations.  During  the  1991/1992  academic  year,  Professor 
Gad-el-Hak  was  a  visiting  professor  at  the  Institut  de  Mecanique  de  Grenoble,  France. 
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RESUME 


Ho,  Chlh-Mlng 

Professor  ho.  Chili-Ming  received  Ills  Ph.D.  In  Mechanics  from  The  Joints  Hopkins  University 
and  his  BS  in  Mechanical  Engineering  from  National  Taiwan  University.  He  is  a  Professor  of  the 
Mechanical.  Aerospace  and  Nuclear  Engineering  Department  at  University  of  California,  Los 
Angeles.  Dr.  Ho  is  the  Director  of  Center  for  Advanced  Technology  of  Micro  .Systems  (('• 
ATOMS).  He  was  a  Professor  of  Aerospace  Engineering  at  the  University  of  Southern  California. 
He  woo  elected  Fellow  of  the  American  Physical  Society  for  pioneering  research  In  controlling  the 
evolutions  of  turbulent  free  shear  layers.  Dr.  Ho  is  known  for  his  work  in  the  areas  of  turbulence, 
aerodynamics,  aeroacoustics  and  bio  fluid  dynamics.  He  also  specializes  in  sensor  development 
and  recently  established  the  C-ATOMS  at  UCLA  for  studying  imcio-elcctto-mcchanicaFsystems. 
Ho  has  published  more  than  eighty  papers  In  these  fields  and  received  US.  patents  in  optical  sensor 
and  in  nozzle  design  to  enhance  mass  transfer.  He  lias  presented  numerous  keynote  talks  as  well  as 
served  as  chairman  in  international  and  national  technical  conferences.  He  is  an  Associate  Editor  of 
the  A5ME  Journal  of  Fluids  Engineering  and  was  an  Associate  F.ditor  of  the  A1AA  Journal  from 
1985  to  1987. 

Dr.  Ho’s  research  projects  include  microchannel  How  with  integrated  microscnsors,  noise 
production  In  turbine  machinery,  mass  transfer  in  turbulent  wakes,  instabilities  in  free  shear  layers, 
flow  in  human  nasal  passage,  supersonic  jet  noise  and  unsteady  aerodynamics. 

HONORS 


Fellow  of  the  American  Physical  Society 

Citation:  For  contributions  in  understanding  the  sensitivity  of  free  shear  layers  under 
perturbations  or  geometrical  v  am  lions  in  initial  conditions.  Ilis  pioneering  works  have 
been  served  as  bases  for  controlling  the  evolutions  of  turbulent  free  shear  layers. 

Honorary  Professor  of  Nanjing  Aeronautical  Institute,  China 

Associate  Fellow  of  the  American  Institute  of  Aeronautics  and  Astronautics 

Phi  Beta  Kappa  Tau  Beta  Pi  Sigma  Xi  Phi  Tau  Phi 

MEMBERSHIP  IN  PROFESSIONAL  .SOCIETIES 

AIAA  ASMR 

ADVISORY  COMMITTEE 

Technical  Committee  of  Aeroacoustics 
American  Institute  of  Aeronautics  and  Astronautics 

Advisor  lu  Institute  of  Aeronautics  and  Astronautics 
National  Cheng* Kung  University,  Tainan,  Taiwan 

Chair  of  Aviation  Technology  Advisory  Committer. 

Center  for  Aviation  and  Space  Technology 
Indue  trial  Technology  Research  Institute,  Hsinr.hu,  Taiwan 

REFERENCES 

Who’s  Who  in  Science  and  Engineering  Who’s  Who  in  Technology  Today 

Who’s  Who  in  Frontier  Science  and  Technology  Who’s  Who  in  the  West 
Who's  Who  of  Emerging  Leaders  in  America 


APS 

1980-84: 

1987-89: 

1992-94: 
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Vitae  of  Yu-Chong  Tai 


Assistant  Professor  of  Electrical  Engineering 
California  Institute  of  Technology 
Pasadena,  CA  91 125 
*  (818)356-8317 


(Pil  Redacted] 


Education 

B.S.  Naiional  Taiwan  University  (1981). 

M.S.  £ECS,  University  of  California,  Berkeley  (1986). 

Ph.D.  EECS,  University  of  California,  Berkeley  (1989). 

Academic  Experience 

1983-1988  Research  Assistant  (EECS),  U.C.  Berkeley. 

1986-1987  Teaching  Assistant  (EECS),  U.C.  Berkeley. 

1988- 1989  Graduate  Instructor  (EECS),  U.C.  Berkeley. 

1989- now  Assistant  Professor  (EE),  California  Institute  of  technology 

Awards  and  Honors 

Ross  N.  Tucker  Award,  AIME  Electronics  Material  Award  (1987). 

Outstanding  Student  Paper  Award,  IEDM  (1988). 

IBM  Fellowship  (1988-1989). 

David  J.  Sakrison  Memorial  Prize,  U.C.  Berkeley  (1989). 

A 

Earnest  C.  Watson  Lecturer,  Caltech  (1990). 

Phi  Tau  Phi  (1990). 

Presidential  Young  Investigaior  (PYI)  Award,  NSF  (1991-1996). 

David  and  Lucile  Packard  Fellowships,  David  and  Lucile  Packard  Foundation  (1991-1996). 
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V.VJA 


Dr.  C.  R.  Smith 

Professor  of  Mechanical  Engineering  and  Mechanics 
Packard  Laboratory,  #  19 
Lehigh  University 
Bethlehem,  PA  18015 
(216)  768-4116 

Education 

B.$.  Mechanical  Engineering,  Stanford  University,  1966 
M.S.  Mechanical  Engineering,  Stanford  University,  1968 
Ph.D.  Mechanical  Engineering,  Stanford  University,  1971 

Professional  Experience 

Assistant  Professor,  Mechanical  Engineering,  Purdue  Univ.,  1971-76 
NAS  A- AS  EE  Faculty  Fellow,  NASA-Ames  Research  Ctr.,  1974 
Associate  Professor,  Mechanical  Engineering,  Purdue  Univ.,  1976-78 
Associate  Professor,  Mechanical  Engineering,  Lehigh  Univ.,  1978-82 
Professor,  Mechanical  Engineering,  Lehigh  Univ.,  1982-present 
Director  of  Industry  Liaison  Program  of  Lehigh  University 
Thermo-Fluids  Institute,  1980-86 

Division  Chairman,  Thermo-Fluids  Division,  ME/Mech,  Lehigh  Univ.,  1982-87 
Research  and  Professional  Activities 


(PH  Redacted] 


•  Has  conducted  research  in  experimental  fluid  dynamics,  turbulent  boundary  layers, 
turbulent  flow  structure,  fluid  mixing,  flame  stabilisation,  flow  visualization,  image 
processing,  vortex-surface  interactions,  wake  flows,  unsteady  separation  phenomena,  and 
"smart”  material-fluid  interaction. 

•  Has  been  principal  investigator  on  contracts  and  grants  totalling  over  $3,600,000. 

•  Has  supervised  seven  Ph.D.  theses  and  seventeen  M.S.  theses  (currently  supervising  three 
Ph.D.  and  two  M.S.  theses) 

•  Co-founder  of  Vortab,  a  company  involved  in  the  development  and  marketing  of 
improved  flow  conditioners  and  static  mixing  devices 

•  Has  consulted  for  over  twenty  firms  (1973-present) 

•  A  member  of  ASME,  Pi  Tau  Sigma,  Tau  Beta  Pi,  American  Physical  Society,  and  AlAA 


Bm graphical  Sketch 

Dr.  Smith  is  a  Professor  of  Mechanical  Engineering  at  Lehigh  University,  Bethlehem, 
Pennsylvania,  and  is  internationally  recognized  for  his  research  in  fluid  mechanics, 
particularly  the  areas  of  turbulent  flow,  unsteady  flow  effects,  flow  visualization  techniques,  fluid 
dynamic  image  processing,  flow  mixing  and  conditioning,  and  fluid  dynamic  use  of  "smart” 
materials.  He  has  published  over  seventy  five  archival  journal  and  proceedings  papers,  and  has 
made  over  one  hunderd  conference,  workshop,  and  invited  seminar  presentations. 


A-ll/A-12 
Reverse  Blank 


Demetri  P.  Telionis 


PERSONAL  DATA  EDUCATION  _ 

(Plj  Redacted] 

Cornell  University  (1970)  Ph.D,  in 
Aerospace  Engineering 
Cornell  University  (1069)  M.S.  in 
Theoretical  &  Applied  Mechanics 
National  Technical  University  of  Greece 
ee-(703l  231-7492  (1964)  Diploma 


MILITARY  SERVICE  Served  as  Ensign  in  the  Royal  Greek  Navy  from  1964  to  1967. 

HONORS  AND  RECOGNITIONS  The  McMullen  k  Glauert  Fellowships  at  Cornell 
University.  Invited  keynote  speaker  at  (i)  CAS1/AIAA  Joint  Aeronautical  Meeting, 

!ii)  AFOSR  International  Symposium  on  Unsteady  Aerodynamics,  1976;  (in)  AGAkD 
nternational  Symposium  Unsteady  Aerodynamics  1977:  f iv )  AGARD-FDr-VKI  Special 
Course  Series,  1980. 

The  College  of  Engineering  teaching  award  1979.  The  SE  American  Association  of 
Engineering  Education  Research  Award  for  best  paper  in  1987.  The  ESM  Frank  E.  Maher 
Outstanding  Educator  Award  1990. 

MEMBERSHIP  QELEARNEP  SOCIETIES 

AIAA  (Associate  Fellow),  ASME,  IUTAM,  ASEE,  Sigma  Xi,  Virginia  Academy  of  Science, 
EDITORIAL  WORK 

Associate  Editor  of  the  ASME  Journal  of  Fluids  Engineering  1987—90 
Editor  of  the  ASME  Journal  of  Fluids  Engineering  1990-95 

PROFESSIONAL  EXPERIENCE 

GREEK  NAVY  SHIPYARDS,  January  1965  —  February  1966,  Mechanical  Engineer 

{while  serving  as  ensign  officer) 

5M,  Binghampton,  NY,  USA,  June  1968  -  October  1968,  Consultant. 

BMT  Boatyard,  Athens,  Greece,  May  1977  -  June  1981,  Consultant 

FUSION  SYSTEMS,  INC.,  Washington,  D.C.,  July  1986  -  Dec.  1989,  Consultant. 

ITT,  INC.,  Roanoke,  VA.,  August  1989  -  August  1990,  Consultant. 

TEACHING  EXPERIENCE 

VIRGINIA  POLYTECHNIC  INSTITUTE  AND  STATE  UNIVERSITY,  Blacksburg,  Va 
September  1970  —  August  1974,  Assistant  Professor 
September  1974  —  August  1978,  Associate  Professor 
September  1978  —  Present,  Professor 

INVITED  TALKS  AND  SEMINARS  18  invited  seminars  presented  at  universities  or 
research  institutes 
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PRESENTATIONS  About  76  presentations  at  national  and  international  meetings. 
PUBLICATIONS 

Book:  pnsteady  Viscous  Flow.  Springer  Verlag,  New  York,  Heidelberg,  Berlin,  Nov.  1981. 

Papers:  Approximately  80  referred  publications  in  Journals  and  Proceedings  (at  least  70  of 
which  are  lull  papers). 

FUNDED  RESEARCH  Directed  and  participated  in  funded  research  in  the  past  10  years 
for  NASA,  NSF,  ONR,  AFOSR.  DOE,  NAVAIR,  ARO,  DOT  (responsible  for  over 
$2,000,000  in  total  expenditures),  in  areas  of  Unsteady  Boundary  Layers— Numerical 
Analysis,  Boundary— Layer  Separation,  Turbulent  Boundary  Layers  and  Closure  Models, 
Homogeneous  Shear  Turbulent  Flos,  Unsteady  Aerodynamics-Unsteady  Stall,  Flow 
Visualization  Methods,  Unsteady  Waves,  Aerodynamics  of  Ground  Transportation, 
Laser-Doppler  Velodmetry,  Three-Dimensional  Separation,  Vortex  Shedding  Over  Single 
and  Multiple  Bodies,  Blade-Vortex  Interaction.] 
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Anthony  Leonard 
Biographical  Sketch 

Dr.  Leonard’s  area  of  interest  Is  computational  fluid  dynamics  and  its  application  to  a  wide  variety 
of  flows  including  turbulence,  transitional  flows,  and  bluff-body  aieodynamlcs.  The  computational  tools 
he  has, been  concerned  with  are  Lagrangian  vortex  methods,  spectral  methods,  and  the  technique  of 
large-eddy  simulation.  He  has  also  been  Involved,  more  recently,  in  the  application  of  dynamical  sys¬ 
tems  theory  to  fluid  transport  and  mixing.  i 1 

He  studied  mechanical  engineering  at  the  California  Institute  of  Technology  and  obtained  a  B.S. 
degree  in  1959,  He  attended  graduate  school  at  Stanford  University,  receiving  an  M.S.  in  Engineering 
Science  in  1960  and  a  Ph.D  in  Mechanical  Engineering  in  1963.  His  research  experience  in  a 
gove mm ent/indus trial  setting  includes  three  years  and  the  RAND  Corporation  and  twelve  years  at 
NASA’s  Ames  Research  Center  in  the  Computational  Fluid  Dynamics  Branch,  Dr.  Leonard’s  academic 
career  includes  seven  years  at  Stanford  as  Assistant  and  Associate  Professor  of  Mechanical  Engineering 
and  one  year  as  Visiting  Professor  of  Aeronautics  at  Caltech.  For  the  past  six  years  he  has  been  Profes¬ 
sor  of  Aeronautics  at  Caltech. 
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DISTRIBUTION  LIST 


Internal 

Codes  01 X  (Chapman) 

0251  (Conforti) 

0261 
0262  (2) 

02244 

10  (Roderick) 

102  (Lima) 

103  (Sullivan) 

1 06  (Dickinson) 

213  (Strawderman,  Schloemer,  Snarski,  Tito,  Payne) 

214  (Bakewell,  Keith,  Hurdis,  Chen,  Patton,  Berliner,  Abraham) 

22 

33A  (Katz) 

2313  (Keefer) 

411  (Kalinowski,  Wagner,  Kieronski) 

432  (Straw) 

441  (Borgeson) 

80  (Dence) 

8141  (Galib) 

82  (White) 

821  (Nadolink) 

8211  (Aaron) 

8212  (Shonting) 

8214  (Brown,  Philips,  Dick,  Koskie) 

8219  (Goodrich) 

823  (Meng) 

8232  (Hendricks,  Temple,  Hrubes,  Henoch) 

8233  (Uhlman,  Kirschner,  Kaufman,  Marshall,  Kuklinski,  Grant,  Rice,  Myers,  Huyer, 

Krol,  Koch,  Varghese) 

8234  (Kemp,  Jilling,  Hrubes,  Henoch,  Lee,  Doyle,  Savoie,  Bandyopadhyay  (30)) 

83  (McNally) 

8321  (Sjoblom,  Barker,  Parker) 

8322  (Lefebvre,  Cipolla,  Jordan,  Schwemin) 

8323  (Bitsakis) 

Total:  104 


